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ABSTRACT 

We have identified 7824 radio sources from the 1.4 GHz NRAO VLA Sky Survey (NVSS) 
with galaxies brighter than K — 12.75 mag. in the Second Incremental Data Release of the 
6dF Galaxy Survey (6dFGSDR2). The resulting sample of redshifts and optical spectra for 
radio sources over an effective sky area of 7076 deg'^ (about 17 per cent of the celestial sphere) 
is the largest of its kind ever obtained. NVSS radio sources associated with galaxies in the 
6dFGS span a redshift range 0.003 < z < 0.3 and have median z — 0.043. Through visual 
examination of 6dF spectra we have identified the dominant mechanism for radio emission from 
each galaxy. 60 per cent are fuelled by star-formation and 40 per cent are fuelled by an active 
galactic nucleus powered by a supermassive black hole. We have accurately determined the 
local radio luminosity function at 1.4 GHz for both classes of radio source and have found it to 
agree well with other recent determinations. From the radio luminosity function of star-forming 
galaxies we derive a local star formation density of 0.022 ± 0.001 Mq yr^^ Mpc^"^, in broad 
agreement with recent determinations at radio and other wavelengths. 

We have split the radio luminosity function of radio-loud AGN into bins of absolute K- 
band magnitude (Mk) and compared this with the underlying X-band galaxy luminosity func- 
tion of all 6dFGS galaxies to determine the bivariate radio- A'-band luminosity function. We 
verify that radio-loud AGN preferentially inhabit the brightest and hence most massive host 
galaxies and show that the fraction of all galaxies which host a radio-loud AGN scales as 



/radio — loud Z/^ for /radio — loud 

hole mass and stellar mass. 
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< 0.3, indicative of a similarly strong scaling with black 



galaxies: luminosity function, mass function - galaxies: starburst - 



1 INTRODUCTION 

In recent years a new generation of radio surveys has been released 
covering a large fraction of the celestial sphere down to flux den- 
sity levels of a few mjy (eg. NRAO VLA Sky Survey (NVSS), 
ICond on et al.' (19981); Sydne y University M olonglo Sky Survey 
(SUMSS), Bock et al] j 19991) fe lMauch et a ll [2003): Faint Images 
of the Radio Sky at Twenty-cm (FIRST), Becker et al. ( 1995); West- 
erbork Northern Sky Survey (WENSS). iRengelinki ([1998)). The vast 
majority of the radio sources in these surveys are produced by Ac- 
tive Galactic Nuclei (AGN) powered by supermassive black holes 
in galaxies with a median redshift of z ~ 0.8 ( Condon et al.. 1998i) . 
These dominate the radio source population above flux densities of 
lOmJy. Below lOmJy the surveys contain an increasing fraction of 
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nearby (z < 0.1) ga laxies whose r adio emission is fuelled by on- 
going star formation ( ICondoiJll989l) . Because radio surveys probe a 
wide range of redshifts, studying their global properties statistically 
(eg. through radio source counts or luminosity functions) provides a 
powerful constraint on the evolutionary properties of massive galax- 
ies throughout the history of the universe. Unfortunately, radio sur- 
vey data alone are not sufficient to constrain the current models of 
radio source counts. These models are strongly dependent on the 
form of the local radio luminosity function (RLF) from which the 
mea sured source counts can be extrapolated via evolut ionary models 
(eg. lDunloD&Pea^ll99d : I Jackson & Wall 1 19991) or compared 
directly to the measured radio luminosity function of samples se- 
lected at higher redshift (eg.^rown et al. 2001; Sadler et al. 200^. 
The current ge neration of r edsh ift surveys (eg. Sloan Digital Sky 
Survey (SPSS), 'York et alj ^20001); 2 d egree Field Galaxy Redshift 
Survey (2dFGRS), CoUes s et aimOOlh ; 6 degree Field Galaxy Sur- 
vey (6dFGS), Jones et all ( 20041) ) provide a useful tool for calculat- 
ing the local radio luminosity function, as they provide redshifts for 
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thousands of radio sources in the local universe. The host galaxy 
spectra provide additional value in this context as they can be used 
to determine the physical cause of the radio emission from galaxies, 
thereby disentangling the star-forming galaxy population from the 
radio-loud AGN. 

Table □ compares the underlying radio-optical samples used 
in recent measure ments of the local radio luminosity function. 
ISadler et alJ l l2002h identified NVSS radio sources in the 2dFGRS 
and derived the local radio luminosity function of both star-forming 
galaxies and radio-loud AGN; this will hereaf t er be r efen'ed to as the 
2dFGRS-NVSS sample. Magliocchetti et al.l ( |2002|) identified opti- 
cal counterparts to FIRST galaxies in the 2dFGRS; this will here- 
after be referred to as the 2dFGRS-FIRST sample. Both the radio- 
selected 2dFGRS samples cover a relatively small area of sky to red- 
shifts of 2: ~ 0.3 and predominantly comprise radio-loud AGN. In 
particular the 2dFGRS-FIRST sample detected relatively few star- 
forming galaxies as the high-resolution FIRST survey resolves out 
much of the ext ended radio emission in the disks of nearby galaxies. 
iBest et alj JioOSh) identified NVSS and FIRST radio sources in the 
second data release of the SDSS; this will hereafter be referred to as 
the SDSS-NVSS/FIRST sample. They optimised their radio-source 
selection via a hybrid method which used both the resolution of the 
FIRST survey and the surface brightness sensitivity of the NVSS 
survey. However, they detected relatively few star-f orming galax- 
ies ab ove their adopted flux density limit of 5 mjy. ICondon et alj 
( |2002|) crossmatched the NVSS catalogue with galaxies brighter 
than mp = 14.5 mag. in the Uppsala Galaxy Catalogue (UGC; 
lNilsorJll973h and constructed a a sample of 1966 radio sources in 
the local [z < 0.03) universe across 4.33 sr of the northern sky; 
this will hereafter be referred to as the UGC-NVSS sample. The 
wide yet shallow UGC-NVSS sample contained a large population 
of star-forming galaxies but relatively few radio-loud AGN. Larger 
numbers of radio-loud AGN are found at higher redshifts. 

This paper presents the local radio luminosity function derived 
from a well-defined subsample of the 7824 radio sources from the 
NVSS catalogue identified with galaxies observed in the Second In- 
cremental Data Release (DR2) of the 6dFGS; this will be referred 
to as the 6dFGS-NVSS sample. Galaxies in the 6dFGS-NVSS sam- 
ple presented in this paper lie at redshifts intermediate between the 
more distant 2dFGRS-NVSS and SDSS-NVSS/FIRST samples and 
the nearby UGC-NVSS sample, in an effective sky area much larger 
than that of the 2dFGRS and SDSS derived samples. 6dFGS-NVSS 
galaxies are contained in a larger volume of space than any previ- 
ous local radio source sample, with significant populations of both 
radio-loud AGN and star-forming galaxies all contained in the one 
survey. 

The radio luminosity function derived in this paper has a num- 
ber of advantages over previous determinations: 

• The Near-Infrared (NIR) and radio input catalogues are from 
homogeneous surveys with a single instrument, and avoids biases 
which may result from combining surveys. 

• The digital Jf-band magnitudes from the 2MASSXSC are 
more accurate than those which have been derived from measure- 
ments of photographic plates (eg. 2dFGRS, UGC). 

• The near-infrared selection of the 6dFGS means it is relatively 
unaffected by dust in both the target galaxy and our own Galaxy. 
This means that the radio selected sample will not be biased with 
respect to the amount of dust in the host galaxy, which can steer 
surveys away from galaxies with higher star-formation rates. Near- 
infrared magnitudes are also more closely linked to the old stel- 
lar population in galaxies, effectively tipping the balance in favour 



of massive early type galaxies which preferentially host radio-loud 
AGN. 

• The large sample volume of the 6dFGS contains more radio 
sources than any previous radio-selected galaxy sample (see ta- 
ble [B. 

The structure of this paper is as follows. Section |2] describes 
the 6dFGS and NVSS surveys and our method of identifying radio 
sources in the 6dFGS. Section |3] outlines the global properties of 
6dFGS-NVSS galaxies. Section |4] examines the population of pri- 
mary sample objects which are detected in the IRAS-FSC, both as a 
consistency check for the spectroscopic classification of the sample 
and to derive the radio-FIR correlation for a larger sample of star- 
forming galaxies than any obtained to date. Section |5] presents the 
local radio luminosity function at 1.4 GHz for the radio-loud AGN 
and star-forming galaxies. Section 6 presents a derivation of the star- 
formation density at the present epoch. Section |7] presents the bi- 
variate radio-NIR luminosity function of radio-loud AGN. Finally, 
Section[8]summarises the main results. Throughout, if not explicitly 
stated, we adopt a ACDM cosmology with parameters Sim ~ 0.3, 
Q.A = 0.7, andi^o = 70kms'^Mpc"^ tSpergel et al...2003il . 



2 SAMPLE SELECTION 

2.1 The 6 degree Field Galaxy Survey 

The 6 degree Field Galaxy Survey is a spectroscopic survey of 
about 170000 objects in the 17 046deg^ of sky with declination 
5 < 0° and galactic latitude |6| > 10°. The survey was carried 
out on the 1.2 metre UK Schmidt Telescope at Siding Spring Obser- 
vatory. The majority of targets (113 988) are near-infrared selected 
to be complete to an apparent magnitude limit of K = 12.75 mag. 
from the Two Micron All Sky Su rvey Extended Source Catalogue 
(2MASS XSC; Ijarrett et all I2OOOI) and comprise the 6dFGS pri- 
mary sample. The remainder of the survey is filled out by a num- 
ber of smaller additi onal target sample s selected in various ways 
at other wavelengths jjones et ai]|2004h . The spectroscopic sample 
considered in this paper consists of the 47 317 objects from the K- 
selected primary sample that have been observed during the period 
2002 January to 2004 Oct ober and released in the Second Incremen- 
tal Data Release (DR2; Jones et alj 12005) of the 6dFGS and that 
overlap with southern part of the NVSS survey (—40° < 5 < 0°). 

6dFGS spectra are obtained through 6arcsec diameter fibre 
buttons which correspond to a projected diameter of 6.8 kpc at 
the median redshift of the survey (z 0.05). 6dF fibres in- 
clude an increasing fraction of total galaxy light for higher red- 
shift galaxies, and therefore galaxies with emission-line nuclei are 
easier to recognise at lower redshift. 6dFGS observations prior to 
2002 October were made with reflection gratings and cover a wave- 
length range of 4000-8400A. Since 2002 October observations were 
made with Volume-Phase transmissive Holographic (VPH) gratings, 
which have provided improved efficiency and data uniformity and 
cover a wavelength range of 3900-7500A. Primary target spec- 
tra typically have resolution R ~ 1000 and signal-to-noise ratio 
(S/N) of 10 pixel" ^ allowing redshifts and host galaxy parameters 
(ie. emission-line ratios) to be readily determined. Redshifts are 
measured f r om th e spectra using the RUNZ package described by 
Ijones et ai]( l2004l) and assigned a quality class Q by human opera- 
tors based on the reliability of the measured redshift. Only Q = 3 
or 4 redshifts are expected to be reliable, Q — 6 redshifts are re- 
served for Galactic sources (see Ijones et al. I [20^ for a complete 
description). 
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Table 1. Comparison of recent radio-selected samples used to derive the local RLF. 



UGC-NVSS 
Condon et al. (2002) 



6dFGS-NVSS 
This paper 



2dFGRS-NVSS 
Sadler et al. (2002) 



2dFGRS-FIRST 
Magliocchetti et al. (2002) 



SDSS-NVSS/FIRST 
Best et al. (2005b) 



No. Galaxies 
Mag. limit 
Slim (mJy) 
Area(Sr.)^ 
z range'' 
median z 
Volume (10SMpc3)<: 
T (Gyr)d 
Mpc^^ 



PSF Mq yr" 



Star-forming galaxies 



1672 
ip < 14.5 
2.5 
4.33 
< 0.04 
0.012 
7 

0.543 
0.018 



4006 
Ktot < 12.75 
2.8 
2.16 

<o.i 

0.035 
53 
1.30 
0.022 



242 

14.0 <bj < 19.4 
2.8 
0.10 
< 0.15 
0.043 

8 
1.89 
0.031 



177 

14.0 <bj < 19.45 
1.0 
0.07 
£0.25 
0.100 
22 
2.94 



497 

14.5 < r < 17.77 
5.0 
~ 0.4<= 
<0.15 
0.055 
~ SO'' 
1.89 



No. Galaxies 
Mag. limit 
Slim (mJy) 
Area(sr.)^ 
2 range'' 
median z 
Volume (lO'^Mpc^)" 
T (Gyr)d 



Radio-loud AGN 



294 
rrip < 14.5 
2.5 
4.33 
< 0.04 
0.019 
7 

0.543 



2661 
Ktot < 12.75 

2.8 

2.16 
<0.2 
a073 

391 

2.44 



420 

14.0 < bj < 19.4 
2.8 
0.10 
< 0.3 
0.140 
58 
3.44 



372 

14.0 <bj < 19.45 
1.0 
0.07 
<0.3 
0.150 
40 
3.44 



2215 
14.5 <r < 17.77 
5.0 
~ 0.4" 
< 0.3 
0.165 
~ 230" 
3.44 



NOTES: 

For the 2dFGRS-NVSS and 6dFGS-NVSS samples the area is the effective area as defined in Section [23] 
'' The maximum redshift of each survey has been estimated from the redshift histogram. 
" Volumes are calculated to the maximum redshift of each survey and reduced by the fraction of sky surveyed. 

T i s the lookback tim e in Gyr and is calculated from the maximum redshift of each survey. 
" As lBest et all J2005tl) did not calculate the sky area of the SDSS-NVSS/FIRST sample, we have roughly estimated it from coverage plots. 



The total X-band magnitudes (Ktot) used for selection 
of 6dFGS primary targets have been derived from isophotal 
magnitudes (ifiso) and diameters to an elliptical isop hote of 
fXK = 20 mag. arcsec^^ (equation 1 of I Jones et al.l 120(3) . These 
derived values of Ktot are more robust than those quoted in the 
2MASS XSC and have therefore been used throughout this paper. 
The derived K magnitudes have typical errors AK < 0.1 mag. and 
the near-infrared- selected sample is complete to the survey limit of 
K = 12.75 mag. j Jones et al.ll2006l) . 



2.2 The NRAO VLA Sky Survey 

The NVSS is a radio imaging survey of the entire sky north of 
S = —40° at 1.4 GHz carried out on the Very Large Array (VLA) 
telescope in its DnC and D configurations. Its principal data prod- 
ucts are a set of 2326 4° x 4° continuum images and a catalogue 
of about 2 X 10® fits of elliptical gaussians to discrete sources 
stronger than S'i.4ghz ~ 2mJy. The images have 6 = 45arcsec 
FWHM resolution with position accuracy < 1" for so urces stronger 
than 15 mJy, increasing to > 7" at the survey limit dCondon et al.l 
Il998h . The resolution and positional accuracy of the NVSS makes 
identification of radio sources with objects on optical survey plates 
straightforward in the majority of cases. There are 580419 NVSS 
catalogue sources which overlap with the 6dFGS; a region bounded 
by 0° > (5 > —40°, |6| > 10°. These sources form the primary 
1.4GHz-selected input catalogue for the 6dFGS-NVSS sample. 



2.3 Sky coverage 

The 6dFGS uses a tiling algorithm wi th variable overlap de pend- 
ing on the underlying galaxy density ( ICamDbeIletZll2004l) . Be- 
cause of_this_and because not all 6dFGS tiles were observed in the 
DR2 j Jones et^l2005t) . calculating the exact area of sky surveyed 
is not straightforward. A method of calculating sky area for surveys 
with gaps in coverage such as the 6dFGS DR2 was described by 
Folkes ( 1999) and was also applied to 2dFGRS-NVSS galaxies by 
I Sadler et al.. 12002). This method estimates the area of sky covered 
by dividing the number of galaxies observed by the mean surface 
density of galaxies in the survey target list. The surface density of 
6dFGS primary targets is 6.7 objects deg^^. 47 317 primary targets 
north of —40° have been observed in the DR2. Therefore the ef- 
fective area is 7 076deg'^, or 17. 15 per cent of the celestial sphere. 
Calculation of the effective area in this way also accounts for the 
5 per cent fibreing incompleteness of the 6dFGS. However, it takes 
no account of the small amount of incompleteness due to failed ob- 
servations (eg. low S/N spectra); this is discussed in more detail in 
Section|5] 

2.4 Crossmatching the 6dFGS primary sample with the NVSS 

When searching for NVSS radio source identifications of the 6dFGS 
primary targets we used a method which aimed to maximise both the 
completeness (ie. all potential radio sources are included) and the 
reliability (ie. all included radio identifications are genuine) of the 
database. To ensure a high completeness and reliability of NVSS- 
2MASS identifications we determine a maximum position offset 
which includes all possible identifications, construct a list of these 
identifications and then verify each of them by visually inspecting 



© 2006 RAS, MNRAS Q00.mt2T] 



4 T. Mauch and E. M. Sadler 



1000 : 



100 : 



10 





- 










i 





30 



10 



20 30 40 

Position offset (arcsec) 



50 



60 



Figure 1. For all 25 096 6dFGS primary targets with a single NVSS cata- 
logue source within 3 arcrain, the open circles linked by a solid line show 
the distribution of position offsets to the radio source. The dashed line shows 
the average distribution of position offsets to NVSS catalogue sources from 
positions in 5 random catalogues with the same size and sky coverage as the 
6dFGS primary sample, also with a single radio source within 3 arcmin. 



overlays of radio contours onto optical images. Although checking 
all identifications by eye is time consuming, in our experience it is 
the best way to minimise the number of false identifications, as at 
all separations some radio identifications selected by position offset 
alone will be chance coincidences. Visual inspection of every radio- 
optical overlay also allows us to correct the catalogued flux densi- 
ties of confused or diffuse sources for which the elliptical gaussian 
model used in the NVSS catalogue does not measure reliable radio 
source parameters. Our visual identification m ethod closely follow s 
that described by both . Condon et al.1 j2002h & ISadler et ^ ( l2002h : 
the reader is referred to those papers for examples and descriptions 
of radio source identifications with bright galaxies in optical images. 

To include all possible radio source identifications, we decided 
to consider all radio sources with NVSS-2MASS position offset less 
than 3 arcmin. This offset was chosen as it is large enough to include 
the majority of possible multi-component radio galaxies while still 
smaller than the average separation of unrelated radio-sources in the 
NVSS catalogue. Of the 74 609 6dFGS primary sample targets north 
of -40° , 34 097 of them had at least one identification in the NVSS 
catalogue within 3 arcmin, of which 21 597 were observed in the 
6dFGS DR2. We chose to examine the 15 716 6dFGS DR2 galaxies 
that have a single NVSS identification within 3 arcmin separately 
from the 5 881 with more than one. 



2.4.1 Single-component radio sources 

For all 25 096 6dFGS primary targets with a single NVSS catalogue 
match within 3 arcmin, the solid line in Fig.[T]shows the distribution 
of position offsets of radio sources from the 6dFGS galaxy (in arc- 
sec). Five random catalogues containing 74 609 positions (the same 
as the number of all 6dFGS primary targets) and covering the same 
region of sky as the 6dFGS-NVSS sample were also matched with 
the NVSS catalogue. For the random positions with only a single 
NVSS match within 3 arcmin, the average of the 5 resulting dis- 
tributions of NVSS-2MASS position offsets is shown as a dashed 
line in Fig.[T] At 40 arcsec the expected number of random matches 
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Figure 2. The results of visual identification of 5 829 6dFGS DR2 primary 
sample objects with a single NVSS catalogue match within 3 arcmin. The top 
figure shows the 5 016 real identifications, we found no genuine identifica- 
tions with NVSS-2MASS offset > 30 arcsec. In the bottom figure, diamonds 
show the 693 false identifications and stars show the 120 uncertain identifi- 
cations. The dashed line shows the adopted cutoff used to accept uncertain 
identifications. 



is the same as the number of real matches so we conclude that all 
single-component NVSS matches with 6dFGS primary targets that 
have position offset > 40 arcsec are not genuine. 

5 829 6dFGSDR2 primary targets with a single NVSS cat- 
alogue match within 3 arcmin have position offset < 40 arcsec. 
The dashed line in Fig. [T] implies that at separations < 40 arcsec 
some misassociations are expected. We visually inspected overlays 
of NVSS ra dio contours onto bj images from the SuperCOSMOS 
Sky Survey faamblv et al]|200lh for each candidate identification 
and classified them as follows 

• Genuine (5 016 candidates), for objects whose radio-optical 
identification is unambiguously real. 

• False (693 candidates), for objects whose radio-optical identi- 
fication is unambiguously not real. This is usually because the radio 
source is identified with another galaxy on the optical image. 
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Figure 3. Radio contours from the NVSS overlaid onto bj images from tlie SuperCOSMOS Sky Survey for two accepted 6dFGS-NVSS identifications with 
fitted positions in the NVSS catalogue more than 15 arcsec from galaxy positions in the 2MASS XSC. In both images crosshairs show the 2MASS XSC galaxy 
position and a small cross labels the fitted position of the radio source in the NVSS catalogue. The contour levels plotted in both images are: 2.5, 3, 4, 5, 7, 9, 12, 
16, 21, 26 & 30mjy beam^^ . The left image shows the 6dFGS galaxy g2 155362- 150424 aX z = 0.07. The catalogued Si, 4 ghz = 55 mjy NVSS radio source 
is 17.5 arcsec from the host galaxy position. Its Aa spectrum (as defined in Table|2j implies that this is a radio-loud AGN and all of the extended radio emission is 
associated with the host galaxy. The right image shows the 6dFGS star-forming galaxy g2259267-374210 at 2 = 0.004. The catalogued 51.4 qhz = 39.8 mJy 
NVSS radio source is 18.9arcsec from the host galaxy position. This source has a complex radio structure and we have therefore corrected its 1.4GHz flux 
density to 42.5 mJy by adding up the pixels in the radio image. 



• Uncertain (120 candidates), for objects whose radio optical 
identification is ambiguous. This is usually because two possible 
host galaxies can be seen within the NVSS position error ellipse. 

Fig. [2] shows our hand classification of 6dFGS-NVSS candi- 
dates. We find no genuine identifications of single component radio 
sources with position offset between 30 arcsec and 40 arcsec. The 
excess of real over random identifications of NVSS radio sources 
at these separations seen in Fig. [T] is likely to have been caus ed by 
the clustering of galaxies as was noted by iBest et alJ ( l2005bl) : at a 
real galaxy position, there is more likely to be a companion which is 
the true radio source identification within 100 arcsec than at a ran- 
dom position on the sky. For accepted NVSS-2MASS identifica- 
tions the increasing range of position offsets with decreasing NVSS 
flux density re flects the increase in position uncertainties for fainter 
radio sources l lCondonetalll 19981) . It is clear that there is no sim- 
ple NVSS-2MASS offset cutoff which separates real identifications 
from false identifications. The vast majority of real identifications 
have small (ie. < 15 arcsec) position offsets though there is still 
about 2 per cent contamination by false identifications with smaller 
separations at all NVSS flux densities. Conversely, some real identi- 
fications are found with offset > 15 arcsec. This is due to a number 
of factors: 

• Complex and asymmetric radio sources can have fitted posi- 
tions in the NVSS catalogue which differ by more than 15 arcsec 
from the galaxy positions in the 2MASSXSC. Two examples are 
shown in Fig. [3] 

• Some galaxies can appear at the position of one of the lobes 
of a double radio source whose centroid is associated with another 
galaxy. This can lead to false identifications of radio sources with 
offsets < 15 arcsec. 

• One of a pair of galaxies with small separation can be identified 



with a radio source which is clearly identified with the other galaxy. 
This can lead to false identifications of radio sources with offsets 
< 15 arcsec. 

We found 120 candidate identifications which had two or more 
optical galaxies visible on SuperCOSMOS bj images inside the 
3(7 NVSS position error ellipse (shown as stars in Fig. [2}- These 
candidates were considered too uncertain to make a definitive vi- 
sual classification. We chose to accept as a genuine match the 55 
objects with position offset < 20 arcsec for radio sources with 
■51.4 GHz < 5 mJy and position offset < 5 arcsec for radio sources 
with 51.4 GHz > 5mJy. This cutoff is shown as a dashed line in 
Fig. |2] In the region of the figure below this line more than 99 per 
cent of unambiguous identifications (ie. both genuine and false) are 
genuine, so we expect that less than 1 per cent of unsure identifica- 
tions will be misidentified by applying the cutoff. 

2.4.2 Multiple-component radio sources 

A further 5 881 6dFGSDR2 galaxies had more than one identifi- 
cation in the NVSS catalogue within 3 arcmin. We chose to exam- 
ine all such candidates by visual examination of SuperCOSMOS bj 
images overlaid with NVSS contours. In 3 226 of them, one of the 
radio components lay within 40 arcsec of the galaxy. These were 
either the cores of multiple-component radio sources or a single 
component radio source and other unrelated NVSS radio sources 
within 3 arcmin. We chose to accept or reject such candidates us- 
ing the method outlined in Section l2.4. ll for single component radio 
sources and accepted further 2716 radio source identifications. We 
determined NVSS flux densities for accepted multiple-component 
radio sources by summing the individual catalogued flux densities 
of each component. 

We visually classified the remaining 2 655 candidate multiple- 
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component radio sources which have no NVSS catalogue match 
within 40 arcsec of the host galaxy. Classification of such objects 
by eye is straightforward as the sensitivity of the NVSS to ex- 
tended radio emission makes it easy to identify large connected 
structures in the radio images. For brighter radio sources, we sup- 
plemented our visual classification by double-checking our identifi- 
cations with previous optical identifications of radio sources from 
the Revised Third Cambridge catalogue of radio sources (3CR; 
Laing et al. 1983) , the P arkes-MIT-NRAO survey catalogue (PMN; 
Griffith &WrigM 1 1993h and the Molong lo Reference Catalogue 
of radio sources (MRC: lLargeetalJll98li) . A further 37 multiple- 
component radio sources were identified out of the 2 655 candidates. 

2.4.3 Giant radio galaxies 

A detection of at least one NVSS radio component within 3 arcmin 
of the host galaxy is required to identify a radio source using our 
identification procedure. This implies that we miss any popula- 
tion of Gi ant Radio Galaxies (GRGs) I Schoenmakers et al]|200ll : 
ISaripalli e t al. 2005) which have linear sizes in excess of 6 arcmin 
in radio images and no detecti on of a core. About half of the 
GRGs in a sample selected by ISaripalli et al. I |2005D have lobes 
that are detected at the limit of the SUMSS catalogue but cores 
that are not. The space density of GRGs with linear size greater 
than 5 arcmin measured by Saripalli et al. (2005) is 10^^ Mpc^'^ for 
z < 0.4, which implies that about 40 such GRGs would be found 
in the 6dFGS-NVSS sample volume. About 1/3 of these would be 
missed by our identification procedure (ie. those with linear size 
> 6 arcmin and no detection of a core). GRGs are powerful radio 
sources and have typical 1.4 GHz radio powers between 10^^'^ and 
lO^^WHz^^. We expect that the incompleteness of the 6dFGS- 
NVSS sample in this radio power range will be no greater than 
1 percent because of missing GRGs. 

2.4.4 Reliability and completeness 

It is possible to estimate the reliability and completeness of the 
database for single component radio sources using the data pre- 
sented in Fig. [T] We expect our visual classification scheme to be 
most robust for identifications of radio sources within 15 arcsec of 
the host galaxy position. Integrating under the dashed line in Fig. [T] 
out to 15 arcsec and scaling to the number of objects observed im- 
plies that 81 random associations with radio sources are expected. 
During visual classification we rejected 74 candidate matches and 
therefore an estimated 7 accepted identifications in this region are 
spurious. The remainder of the accepted matches had offsets be- 
tween 15 arcsec and 30 arcsec. We expect 238 random associations 
with offsets in this region and during visual identification rejected 
361 candidate matches. This implies that 123 genuine identifications 
in this region have been spuriously classified as false. We believe 
that the majority of these spuriously rejected identifications have 
been made for fainter (Si.4ghz < 2.8 mJy) radio sources many 
of whose NVSS position errors exceed 10 arcsec. Finally, we ac- 
cepted no radio source identifications with position offset between 
30 & 40 arcsec, however based on the small excess of real over 
random matches seen in this region of Fig. [T] we expect there to 
be 51 genuine associations. This means that for single-component 
radio sources in the final database, 7/5016 accepted matches are 
spurious, corresponding to a reliability of close to 100 per cent, and 
(123+51)=174 matches have been erroneously rejected, correspond- 
ing to a completeness of 97 per cent. Our estimate of completeness 



Table 2. Spectral classes visually assigned to 6dFGS-NVSS objects. 



Class 


Type of spectrum 


Number 


Aa 


Pure absorption line spectrum. 


2124 


Aae 


NaiTow LINER-like emission lines. 


370 


Ae 


Emission line AGN. 


370 


SF 


Star-forming galaxy. 


4625 


7 


Unclassifiable or unknown. 


321 


Star 


Galactic star. 


14 


Total 




7824 



is likely to be a lower limit because many of the excess real over 
random matches at separations greater than 15 arcsec are associ- 
ated with galaxy pairs as was noted in section l2.4. II The reliability 
and completeness of the 2716 identifications with more than one 
NVSS catalogue identification within 3 arcmin and a radio compo- 
nent within 40 arcsec of the host galaxy mimics that of the single- 
component radio sources. 

The reliability and completeness of multiple-component radio 
source identifications with no NVSS detection with 40 arcsec of 
the host galaxy is difficult to determine. For the brightest sources, 
we crossmatched the 6dFGS-NVSS database with well-studied 
databases of optical identifications from the 3CR, PMN and MRC 
databases and found no spurious or missing identifications which 
implies our completeness and reliability are 100 per cent. We as- 
sume this extends to fainter radio source identifications in the 
6dFGS-NVSS sample. Multiple-component radio source identifica- 
tions are rare (they comprise < 1 per cent of all radio sources in this 
sample) and we believe these cannot affect the completeness and 
reliability in a significant way. The missing population of GRGs de- 
scribed in Section [2.4.3| also contributes < 1 percent to the incom- 
pleteness of 6dFGS-NVSS identifications. 

In summary, we find that the completeness of the 6dFGS- 
NVSS database is better than 96 per cent (a 3 per cent contribution 
from missing single-component matches and a 1 per cent con- 
tribution from missing multiple-component matches). T he reliabil- 
ity of t he database is better than 99 per cent. Recently iBest et al] 
( l2005bl) described a robust identification method utilising both 
the FIRST and NVSS surveys to find radio sources in the SDSS 
and their matching procedure had an estimated completeness of 
94.4 per cent and reliability of 98.9 per cent. The FIRST survey only 
covers about 1 per cent of the region observed in the 6dFGS DR2 
and we were therefore unable to use their method. We achieve sim- 
ilar (or slightly better) completeness and reliability to them by vi- 
sually classifying all radio-optical identifications; we believe this is 
the most robust method of identifying radio sources in optical sky 
surveys when only the NVSS is available. 

2.5 Classification of 6dFGS optical spectra 

The 6dF spectrum of each accepted radio source identification was 
inspected visually to determine the dominant physical process re- 
sponsible for the rad io emission of each galaxy. We followed the 
scheme described bv lSadleretalJ 1 I2OO2I) and classified each spec- 
trum as either 'AGN' for spectra indicative of galaxies harbouring 
a radio-loud active galactic nucleus or 'SF' for spectra indicative of 
galaxies with ongoing star formation. AGNs can have a pure absorp- 
tion line spectrum like that of a giant elliptical galaxy (classed as 
Aa); a spectrum with absorption lines and weak narrow LINER-like 
emission lines (classed as Aae); a conventional Type II AGN spec- 
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Figure 4. Example 6dFGS spectra of 4 spectral classes. In each spectrum common redshifted absorption/emission lines are labelled at the measured redshift. 
Top left: An 'SF' galaxy, with characteristic strong and narrow Balmer lines. Top right: An 'Aa' galaxy, containing broad Mg absorption as well as a distinct 
break close to the H & K lines. Bottom Left: An 'Aae' galaxy. This looks similar to an 'Aa' galaxy but has a weak [Nil] line. Bottom right: An 'Ae' galaxy. This 
has broad Ho emission and strong [OIII] Unes. 



trum which has nebular emission lines such as [On] , [Olll] or [Nil] 
which are stronger than any hydrogen Balmer emission lines (Ha or 
H/3) (classed as Ae); or a conventional Type I AGN spectrum with 
strong and broad hydrogen Balmer emission lines (also classed as 
Ae). Star-forming galaxies have spectra typical of Hll regions with 
strong narrow emission lines of Ha and H/3 dominating the spec- 
trum (classed as SF). Figure |4] shows an example of each type of 
spectrum. Spectra of galactic stars were classified 'star' and spectra 
for which a classification could not be made (mostly due to a low 
S/N spectrum) were classed '?'. A '?' was also appended to the end 
of any classification which was not certain. Our visual classification 
scheme is summarised in Table|2l 

To check the reliability of the visual classification scheme a 
subsample of ~ 1800 spectra was classified by both of us and the 
classifications compared. The two human classifiers agreed well, 
with disagreement in about 10 per cent of cases. The primary cause 
of disagreement is between the classification of emission-line AGN 
(Aae and Ae classes) and the SF class, and occurs for spectra which 
have line ratios which are borderline between the AGN and SF 
classes. If only the AGN and SF classes are considered then the 
two human classifiers disagree only 5 per cent of the time. In their 
analysis of 2 dFGS spectra using the same classification scheme de- 
scribed here, ISadleretalJfcooa) found that this scheme agre ed well 
with r nethods such as princip al component analysis (PCA; iFolkej 
1 1 9991 ; iMadgwick et al. 2002) and classification based upon diag- 
nostic emission line ratios ( lackson & Londish 2000 ) both of which 
have a similar reliability to that of the visual classification presented 
here. 



3 SAIVIPLE PROPERTIES 



We matched a total of 7 824 6dFGS DR2 primary targets with radio 
sources from the NVSS catalogue, corresponding to a radio detec- 
tion rate of 16.5 per cent. Of the 7 824 detections, 7 672 were iden- 
tified with a single NVSS component and 152 were identified with 
more than one. Table |2] summarises the spectral properties of the 
accepted 6dFGS-NVSS matches. Out of the 7 824 spectra we ex- 
amined, 321 were unclassifiable (classed as '?') and 14 of them 
were Galactic stars (assigned redshift quality Q = 6 in the 6dFGS). 
Of the unclassifiable spectra, 302 had redshift quality Q = 1 or 2, 
these spectra were too noisy for reliable spectral classification. The 
remaining 19 unclassifiable spectra had reliable redshift measure- 
ments, but did not present enough information for reliable visual 
classification. 



The ratio of star-forming galaxies (SF) to radio-loud AGN in 
the sample is approximately 60 per cent SF to 40 per cent radio-loud 
AGN. The fractions of SF to radio-loud AGN in the similarly clas- 
sified 2dFGRS-NVSS sample (Sadler et al. 2002) are 40 per cent SF 
to 60 per cent radio-loud AGN. The difference in these relative frac- 
tions can be explained by the fainter magnitude limit (bj = 19.4) of 
the 2dFGRS, resulting in a larger proportion of detections of more 
distant radio-loud AGN. It is interesting that in both the 6dFGS- 
NVSS and 2dFGRS-NVSS samples, over 70 per cent of all radio- 
loud AGN have absorption-line spectra and would be missed from 
AGN samples selected on the basis of optical emission lines. 
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Table 3. The first 30 entries of the 6dFGS-NVSS data table. 



(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


6dFGS target 


RA 




Dec. 


K 


z 


Q 


Offset 


•S'l.4 GHz 


S'843 MHz 




Spectrum 


name 




(J2000) 












NVSS 


SUMSS 


IRAS FSC 


Class 




{h : m : 


s) (° 


■■'■■") 


(mag.) 






(arcsec) 


(mJy) 


(mJy) 


(Jy) 





g0000124- 


-363113 


00: 


:00: 


:12, 


.39 


-36: 


:31 


:13, 


.1 


12 


.724 


0.1169 


gOOOOMl- 


-251113 


00: 


:00: 


:14, 


.12 


-25: 


:11 


:12, 


.9 


11 


.886 


0.0852 


g0000356- 


-014547 


00: 


:00: 


:35, 


.64 


-01 


:45: 


:47, 


.4 


11 


.267 


0.0246 


g0000523- 


-355037 


00: 


:00: 


:52, 


.33 


-35: 


:50: 


:37, 


.2 


11, 


.330 


0.0521 


g0001197- 


- 140423 


00: 


:01 


:19, 


.74 


-14: 


:04: 


:23, 


.2 


12 


.289 


0.0867 


g0001394- 


-025852 


00: 


:01: 


:39. 


.36 


-02: 


:58: 


:52, 


.1 


12 


.001 


0.1013 


g0001453- 


-042049 


00: 


:01 


:45, 


.29 


-04: 


:20: 


:49, 


.0 


11 


.845 


0.0481 


g0001496- 


-094138 


00: 


:01 


:49, 


.58 


-09: 


:41 


:37, 


.8 


12 


.304 


0.1038 


g0001558- 


-273738 


00: 


:01 


:55, 


.82 


-27: 


:37: 


:38, 


.0 


10 


.164 


0.0283 


g0001567- 


-035755 


00: 


:01 


:56, 


.70 


-03: 


:57: 


:54, 


.8 


12 


.492 


0.0226 


g0001572- 


-383857 


00: 


:01 


:57, 


.20 


-38: 


:38: 


:56, 


.4 


12 


.644 


0.0547 


g0002039- 


-332802 


00: 


:02: 


:03, 


.88 


-33: 


:28: 


:02, 


.2 


11, 


.143 


0.0289 


g0002348- 


-034239 


00: 


:02: 


:34, 


.81 


-03: 


:42: 


:38, 


.6 


10 


.623 


0.0215 


g0002487- 


-033622 


00: 


:02: 


:48, 


.65 


-03: 


:36: 


:21, 


.7 


11 


.455 


0.0208 


g0002547- 


-354319 


00: 


:02: 


:54, 


.69 


-35: 


:43: 


:19. 


.4 


11 


.677 


0.0489 


g0002545- 


-341408 


00: 


:02: 


:54, 


.47 


-34: 


:14: 


:08, 


.4 


10 


.699 


0.0226 


g0002558- 


-265451 


00: 


:02: 


:55, 


.81 


-26: 


:54: 


:51, 


.2 


11 


.994 


0.0665 


g0003051- 


-073700 


00: 


:03: 


:05, 


.06 


-07: 


:37: 


:00, 


.3 


12 


.291 


0.0299 


g0003057- 


-015450 


00: 


:03: 


:05, 


.66 


-01 


:54: 


:49, 


.7 


10, 


.156 


0.0244 


g0003056- 


-295159 


00: 


:03: 


:05, 


.62 


-29: 


:51 


:59, 


.3 


11 


.741 


0.0609 


g0003130- 


-355614 


00: 


:03: 


:12, 


.97 


-35: 


:56: 


:13, 


.4 


10 


.508 


0.0499 


g0003321- 


-104441 


00: 


:03: 


:32, 


.13 


-10: 


:44: 


:40, 


.6 


10 


.348 


0.0299 


g0003449- 


-204757 


00: 


:03: 


:44, 


.90 


-20: 


:47: 


:56, 


.5 


12 


.314 


0.0970 


g0004029- 


-330202 


00: 


:04: 


:02, 


.86 


-33: 


:02: 


:02. 


.0 


12 


.397 


0.0380 


g0004472- 


-013413 


00: 


:04: 


:47, 


.22 


-01 


:34: 


:12, 


.8 


11 


.831 


0.0239 


g0004517-^ 


-060058 


00: 


:04: 


:51, 


.74 


-06: 


:00: 


:57, 


.6 


12 


.722 


0.1080 


g0004576- 


-014108 


00: 


:04: 


:57, 


.57 


-01 


:41 


:07, 


.9 


12 


.569 


0.0239 


g0005026- 


-160715 


00: 


:05: 


:02, 


.58 


-16: 


:07: 


:15, 


.0 


11 


.997 


0.0339 


g0005028- 


-274253 


00: 


:05: 


:02, 


.78 


-27: 


:42: 


:52, 


.5 


11 


.619 


0.0333 


g0005054- 


-070536 


00: 


:05: 


:05, 


.37 


-07: 


:05: 


:36, 


.3 


10 


.065 


0.0128 



4 


9.5 


12.1 




Aa 


4 


3.8 


28.4 




Aa 


4 


15.4 


2.8 




Aae? 


4 


0.0 


48.4 




Aa 


3 


11.0 


40.5 




Aae? 


4 


2.3 


5.5 




Aa 


4 


5.0 


2.8 




Aa 


4 


13.8 


3.4 




Aa 


4 


4.3 


29.8 


0.488 


SF 


4 


2.3 


6.2 


0.466 


SF 


4 


2.5 


11.8 


22.9 


SF 


4 


4.2 


6.2 


0.941 


SF 


4 


0.3 


13.1 


1.121 


SF 


4 


6.9 


3.9 




SF 


4 


24.0 


4.2 




Aa 


4 


2.6 


7.3 


0.760 


SF 


4 


1.1 


94.4 




Aa 


4 


6.3 


3.7 




Aae 


4 


1.1 


7.7 


0.805 


SF 


4 


6.7 


27.8 




Aa 


4 


1.4 


589.5 


1123.4 


Aa 


4 


8.0 


2.5 




Aae? 


4 


18.8 


3.2 


0.274 


SF 


4 


12.9 


2.4 


0.447 


SF? 


4 


5.1 


9.9 


0.843 


SF 


4 


1.5 


9.3 




Aa 


4 


9.3 


3.6 




SF 


4 


13.8 


2.2 


0.268 


SF 


4 


2.9 


10.2 


1.178 


SF 


4 


2.6 


12.9 


1.038 


SF 



Column descriptions: 

(1) The target name of the object from the 6dFGS database. 

(2) & (3) J2000 Right Ascension and declination of the object from the 2MASS XSC. 

(4) Total K-band magnitude calculated from 2MASS isophotal K magnitude accordin g to equation 1 of ljones et al.l j2004h . 

(5) & (6) 6dFGS measured redshift and quality flag as described in jjones et aljj200 4). Only redshifts with Q > 3 were deemed reliable. 

(7) The offset in arcsec from the NVSS radio position to the objects position in the 6dFGS database. 

(8) The 1 .4 GHz flux density in mJy from the NVSS catalogue. 

(9) The 843 MHz flux density in mJy for 6dFGS objects with a SUMSS catalogue match within 10 arcsec. 

(10) The IRAS FSC 60 flux density in Jy for 6dFGS objects which also appear in additional target sample 126 (see ljones et alj2004l) of the 6dFGS. 

(11) The classification of the spectrum as defined in Table[2] 



3.1 The data table 

Table[3]shows details of the sample of 6dFGS-NVSS radio sources. 
The 6dFGS target name is given for each source so that the ob- 
served 6dFGS spectra can be obtained or each obj ect can be matche d 
against other 6dFGS additional target samples jjones et al 
The table also shows 2MASSXSC J2000 positions, derived Tftot, 
redshifts and redshift quahty Q for each source. We list offsets 
from the given position and integrated 1.4 GHz flux densities for 
each NVSS radio source. For single component NVSS radio sources 
with an 843 MHz SUMSS catalogue (version 1.6) detection within 
10 arcsec we also list the 843 MHz flux density of the source. The 
list of SUMSS matches is highly incomplete at present as the 6dFGS 
DR2 covers little of the sky south of 5 = —40° and version 1.6 
of the SUMSS catalogue is not complete north of this declination. 
We plan to thoroughly investigate the properties of SUMSS radio 



sources in the 6dFGS when the full data releases of both surveys 
are made available and this will be the subject of a future paper. For 
objects with a match in the 6dFGS additional target sample selected 
from the Infrar ed Astronomical S atellite (IRAS) Faint Source Cat- 
alogue (FSC) (Mos hir et alJI 19931) (programme id 126; I Jones et al.l 
[2004) we list the 60 pm flux density of the source. The final col- 
umn of the table shows the visually assigned spectral classification 
described in Section [23] 

3.2 Radio source counts 

The differential source counts of 6dFGS-NVSS galaxies are plot- 
ted in Figure |5] These are plotted in a form in which each bin is 
weighted by 5'^''^ such that counts in a static Euclidean universe 
would lie on a horizontal line. Counts determined for the radio-loud 
AGN are plotted as filled circles and those for star-forming galaxies 
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Figure 5. Differential source counts of 1 .4 GHz NVSS sources multiplied by 

5 

S 2 such that counts in a Euclidean universe would lie on a straight Hne. The 
source counts of all 580419 NVSS sources in the region surveyed are shown 
by open squares linked by a solid line. AGN and SF class galaxies in the 
6dFGS-NVSS sample are shown as filled and open circles respectively. The 
eiTor bars shown are the counting uncertainties in each bin. A vertical 
line is plotted at 51.4 gHz = 2.8 mJy which is the completeness limit of 
the 6dFGS-NVSS sample and is also the upper bound of the faintest bins 
calculated from the AGN and SF samples. 



are plotted as open circles. The source counts of all 580419 NVSS 
catalogue sources in the same region of sky as the 6dFGS-NVSS 
sample are shown as squares linked by a solid line. The counts fol- 
low an approximate power law between IJy and 2.8 mJy, below 
which they fall sharply; this is because the NVSS catalogue becomes 
increasingly incomplete below 2.8 mJy ( Condon et al. i998). The 
vertical line at 2.8 mJy shown in the figure divides the two smallest 
bins of the AGN and SF samples. The same falloff can be seen to 
the left of the vertical line in the counts of the AGN and SF subsam- 
ples indicating that this 6dFGS-NVSS sample is complete to a flux 
density limit of 2.8 mJy. 

Radio source catalogues such as the NVSS are known to pre- 
dominantly contain galaxies powered by radio-loud AGN with a 
median redshift of z 0.8 (Condon et al. 1998). In the much more 
nearby (z ~ 0.05) 6dFGS-NVSS sample the source counts of radio- 
loud AGN follow roughly the same power-law slope as that of all 
radio sources but only account for roughly 1 per cent of their total 
number. The source counts of star-forming galaxies stay flat over 
the range 2.8 mJy to 1 Jy and therefore account for only 0.1 percent 
of all sources at 1 Jy rising to 1-2 per cent of all sources at the sur- 
vey limit. The counts of star-forming galaxies and radio-loud AGN 
cross over at about 10 mJy, below which star-forming galaxies dom- 
inate the population of radio sources. At 1 Jy there are over ten times 
more radio-loud AGN than star-formers, but the balance changes at 
lower flux density until at 2.8 mJy there are nearly ten times more 
star-formers than radio-loud AGN in this volume-limited sample. 
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Figure 6. The redshift distribution of the 6dFGS-NVSS sample in bins of 
width 0.005 in redshift. Bottom panel: Redshift distribution of all 7 824 ob- 
jects for which a redshift has been measured. Middle panel: Redshift distri- 
bution of the 4 625 Star-forming galaxies. Top panel: Redshift distribution 
of the 2 864 radio-loud AGN. In each panel a vertical dashed line is drawn at 
the median redshift. 
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Figure 7. The NVSS detection rate of 6dFGS DR2 primary targets. Here the 
detection rate is defined as the percentage of radio identifications of 6dFGS 
DR2 objects for that bin. The error bars shown correspond to the y/n count- 
ing error in each bin. A vertical dotted line is plotted at the median redshift 
of the 6dFGS-NVSS sample and a horizontal dashed line is plotted at the 
overall radio detection rate of 16.5 percent. 



3.3 Redshift distribution 

Fig-Ejshows the redshift distribution of 6dFGS-NVSS objects. The 
redshift distribution of all objects is shown in the bottom panel 
and separate distributions for star-forming galaxies and radio-loud 
AGN are shown in the middle and top panels respectively. The 
median redshift of all radio sources in the 6dFGS-NVSS sample 
is z = 0.046, very close to the median redshift (5 = 0.054) of 
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the underlying if-selected sample jjones et alj|2004h . Star-forming 
galaxies dominate the lower redshift radio source population and 
are found at a median redshift of z = 0.035 while radio-loud AGN 
dominate the higher redshift radio source population at a median 
redshift of 2 = 0.073. 

The optical magnitude limits, radio flux density limits and me- 
dian redshifts of other radio-optical samples, selected in a simi- 
lar way to the 6dFGS-NVSS sample, are presented in Table [T] In 
their analysi s of 6*1 .4 ghz > 2.8 mjy 2dFGRS-NVSS radio sources 
ISadleret alJ ( Il999h found a median redshift of 5 = 0.1 for all 
2dFGRS-NVSS galaxies, z = 0.05 for star-forming galaxies and 
z = 0.14 for radio-loud AGN. The 2dFGRS goes more than two 
magnitudes deeper in bj than the 6dFGS yet the median redshift 
of 2dFGRS-NVSS star-forming galaxies is only slightly higher than 
that found here for 6dFGS-NVSS galaxies, since it is primarily the 
radio flux density limit of optical-radio samples which limits the 
maximum distance to which star-forming galaxies can be found. In 
a survey of 2dFGRS-FIRS T radio sources with a faint er flux density 
limit of Si. 4 GHz > 1 mJv. lMagliocchetti et alj j2002i) found a me- 
dian redshift for star-forming galaxies of 2 = 0.1, more than twice 
that of the NVSS selected samples. Conversely, for radio-loud AGN 
the median redshifts of the 2dFGRS-NVSS and 2dFGRS-FIRST 
samples are roughly similar (z — 0.15) whereas the median red- 
shift of the 6dFGS-NVSS sample is significantly lower. This is a 
consequence of the different optical magnitude limits of the 6dFGS 
and the 2dFGRS; radio-loud AGN fall out of both samples at the 
optical/near-infrared limit of the spectroscopy. 



3.4 Detection rates 

Fig. |7] shows the variation in the detection rate of 6dFGS-NVSS 
radio sources as a function of redshift. The detection rates were cal- 
culated by dividing the number of 6dFGS-NVSS galaxies in redshift 
bins of width 0.01 by the number of the 6dFGS DR2 primary targets 
with Q = 3 or 4 in each redshift bin. The average NVSS detection 
rate is shown as a dashed horizontal line in and the median redshift 
of the 6dFGS-NVSS sample is shown as a dotted vertical line. The 
NVSS detection rate drops sharply from 34 per cent al z — 0.005 to 
just over 11 percent al z — 0.05 (around the median redshift of the 
survey). It then rises steadily to over 20 per cent at z = 0.15 beyond 
which results are dominated by counting errors. This fall and sub- 
sequent rise in the detection rate with redshift can be explained in 
terms of the different redshift distributions of the star-forming and 
radio-loud AGN populations. 

Below z — 0.05, star-forming galaxies dominate the radio 
source population and the detection rate is falling because of the 
1.4 GHz flux density limit of the NVSS. Star-forming galaxies have 
weaker radio powers and fewer are found with increasing redshift 
because the limiting radio power of the 6dFGS-NVSS sample is in- 
creasing with redshift. Above z — 0.05, the detection rate is in- 
creasing and radio-loud AGN dominate the radio source population. 
The population of radio-loud AGN has a higher average radio power 
so they tend not to fall out of the sample because of the flux den- 
sity limit of the NVSS. The K < 12.75 mag. cutoff of the 6dFGS 
means that as redshift increases the galaxies left in the survey are on 
average more luminous in the near-infrared. Radio-loud AGN pref- 
erentially reside within these most luminous galaxies (eg. Fig. [§}, 
resulting in an increase in the detection rate of radio-loud AGN. 
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Figure 8. Bottom: 1.4 GHz radio power vs. Alfc for the 4625 star-forming 
galaxies in the 6dFGS-NVSS sample. Top: 1.4 GHz radio power vs. M/f for 
the 2864 radio-loud AGN in the 6dFGS-NVSS sample. Aa class AGN are 
plotted as dots, Aae class AGN are plotted as circles and Ae class AGN 
are plotted as diamonds. A dashed vertical line is drawn at the value of 
= —24.6 mag. from a Schechter function fit to the local if-band lumi- 
nosity function of galaxies in the 6dFGS by Jones et al. 1 2006). The grey dot- 
dashed lines in each panel show the limiting values of Pi 4 and Mk derived 
from the radio (Si. 4 gHz = 2.8 mJy) and near-infrared (K = 12.75 mag.) 
completeness limits of the 6dFGS-NVSS sample at the labelled redshifts (ie. 
z = 0.015, z = 0.035, the median redshift of star-forming galaxies & 
z = 0.073, the median redshift of radio-loud AGN). These lines show that 
the apparent correlation between Pi. 4 and Mfc is caused by the flux density 
and apparent magnitude limits and the range in redshift of the sample. 

3.5 Luminosity distribution 

Fig. [8] shows the 1.4 GHz monochromatic radio power (Pi. 4, in 
units of WHz^^) vs. the absolute K-band magnitude (Mk) for the 
radio-loud AGN (top) and star-forming galaxies {bottom). Radio 
power measurements have been fc-corrected using the usual form, 
fcradio(z) = (1+2) '^^°' where Q is the spectral iudex (S^ oc v"). 
We determined median spectral indices separately for the AGN and 
SF classes by comparing the flux densities of the 436 galaxies that 
were also found in version 1.6 of the 843 MHz SUMSS catalogue 
jMauch et al]|2003h . For the 228 radio-loud AGN the median spec- 
tral index was fi^ti" = -0.54 ± 0.07 and for the 208 SF galaxies 
the median spectral index was Sgla^ = —0.90 ± 0.07. These val- 
ues of a have been applied to the fe-corrections for all calculations 
of radio power throughout this paper. Both these values are close to 
the value of a = —0.7 often assumed for radio source samples (eg. 
ISadler et alj[2003 : ICondon et al]|2002l) . For most of the galaxies in 
this sample fcradio(z) is negligibly small. 

In the A'-band, fc-corrections {kK{z)) are similar for galax- 
ies of all Hubble types because different amounts of star formation 
in different galaxy types only significantly affect galaxy spectra at 
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wave lengths shorter than 1 /im dCole et al.ll200ll : lGlazebrook et all 
is'-band fc-corrections have been calculated following 
Glazebrook et alj (T99?) who derived fc jc (z) f rom the evolutionary 
synthesis models of Bruzual & Chariot l ll993h assuming an instan- 
taneous burst of star formation at age 5 Gyr. At the small redshifts 
probed by this dataset [z < 0.2), kn^z) never changes Mk by 
more than 0.4 magnitudes. 

The star-forming galaxies in Fig. [8] span a range of Mk be- 
tween —20 and —26, they have median Mk ~ —24.13, lower than 
Mk- Star-forming galaxies, with median log[Pi.4 (WHz^^)] = 
22.13, have weaker radio powers (Pi. 4 < 10^^ W Hz" 1) than radio- 
loud AGN. Radio-loud AGN are almost all found in objects brighter 
than Mk, indicative of their preferential location in the brightest 
galaxies. The median Mk of radio-loud AGN is —25.40 and they 
have median radio power log[Pi.4 (WHz~^)] = 23.04, almost an 
order of magnitude brighter than the median radio power of star- 
forming galaxies. Though star-forming galaxies have lower radio 
powers, radio-loud AGN span a wide range in radio power from 
10^^ to 10^^ WHz-\ indicating that there is no clear observational 
regime in which star-forming galaxies can be separated from radio- 
loud AGN purely on the basis of radio power. 

4 THE FIR-RADIO CORRELATION 

It has been well established that a correlation exists between far- 
infr ared (FIR) and radio continuum emission from normal galaxies 
(seelHelou et al.'(1985[); rCondon & Broderickl (ll988VCo ndon eTal] 
( 11991) : Condon ( 19921)). This correlation has been attributed to on- 
going star formation within the host galaxy, where radio contin- 
uum emission is produced by a combination of synchrotron emis- 
sion from electrons accelerated in the supernova remnants of short- 
lived massive stars (M > 8A/q with lifetimes < 3 x 10*^ yr) and 
free-f ree emission fr om Hll regions ionised by these same massive 
stars ( ICondon|[l992l) . FIR emission is caused by thermal reradia- 
tion of dust in Hll regions heated by this same population of mas- 
sive stars. Though it is understood that recent star formation is the 
process which drives both the radio continuum and FIR emission 
of galaxies which lie on the correlation, the actual mechanism that 
relates non-thermally dominated radio power and thermally domi- 
nated FIR luminosity is poorly understood. 

The FIR-radio correlation has served as a diagnostic tool for 
large samples of radio sources, as it can be used to distinguish 
between galaxies with ongoing star formation and those harbour- 
ing a radio-loud AGN. Galaxies falling on the correlation derive 
their radio emission from star formation, whereas galaxies with 
a radio "excess" above that expected from the correlation derive 
their radio emission from the presence of a rad io-loud AGN (eg. 
ICondon et all |2002| : ICondon & Broderickl Il988h . The correlation 
also constrains models relating star-formation rate to radio emission 
( ICondonl ll992'). In this section results from analysis of the popula- 
tion of radio sources in the 6dFGS-NVSS sample with FIR detec- 
tions in the IRAS FSC are presented, both as a consistency check for 
the spectral classification of the sample and to examine the correla- 
tion for a large sample of normal galaxies. 

4.1 Finding IRAS Detections 

In addition to the primary A'-selected sample, the 6dFGS aims to 
measure redshifts for several additional target samples selected at 
various wavelengths. One of these additional target lists consists of 
~ 11 000 IRAS FSC galaxies selected at 60 from the IRAS FSC 



Table 4. Spectral classification of 6dFGS-NVSS objects which are also 
found in the 6dFGS-FSC sample. 



Spectral class 


,76dPGS-FSC 1 
^'6dFGS-NVSS 


per cent of total^ 


SF 


2690 


58.1 


Ae 


138 


37.3 


Aae 


66 


17.8 


Aa 


4 


0.2 


stai' 


2 


14.3 


? 


42 


13.1 


Total 


2942 


37.6 



NOTES: 

^ The number of 6dFGS DR2 objects which are common to both the 
IRAS FSC and the NVSS catalogue. 

^ The percentage of 6dFGS-NVSS objects that are also 6dFGS-FSC objects. 



(For a complete descript ion of the 6dFGS additional target sam- 
ples see ljones et alj2004l) . Much of this IRAS FSC additional target 
sample (which we call the 6dFGS-FSC sample) overlaps with our 
6dFGS-NVSS sample. It is therefore straightforward to find IRAS- 
FSC detections of 6dFGS-NVSS primary target objects by cross- 
matching their 6dF target ID names. 4403 6dFGS-FSC targets were 
observed in the DR2. Table|4]shows the spectral classification of the 
2 942 galaxies which are common to the 6dFGS-NVSS sample and 
the 6dFGS-FSC sample. Predominantly emission line (ie. Ae, Aae, 
SF) galaxies are detected in the far-infrared as has been observed 
in other sp ectroscopic studies of far- infrared se lected galaxie s (eg. 
iLawrencee t al. 1986: lde Griip et al.| [l992: Sadl er et al|[2002l) . The 
four Aa galaxies detected in the 6dFGS-NVSS-FSC subsample have 
low S/N and may have be en misclassified. 

ICondon et alj ilOoj) noted that the FIR-radio correlation en- 
sures that most far-infrared sources powered by ongoing star for- 
mation are also radio sources and vice versa. The FIR/radio ra- 
tio u — log(5'6o^im/S'i.4GHz) for nearby spiral galaxies in the 
UGC has a mean value of (u) = 2.0 5 ± 0.02 with rms width 
cr„ = 0.2 dCondon & Broderickl [19881) . This value of (u) hap- 
pens to be the ratio between the 280 miy flux density limit of 
the IRAS-FSC at 60 ^im and the limit of the NVSS at 1.4 GHz 
(u ~ 2.05). 58.1 percent of spectroscopically classified SF galaxies 
have 6dFGS-FSC detections leaving 929 SF galaxies with no FIR 
detection. These missing detections are the result of incompleteness 
at the flux density limits of the 6dFGS-FSC and 6dFGS-NVSS sam- 
ples and incom plete sky coverage of the region —40° < 5 < 0° in 
the IRAS FSC ilVIoshir et al.lll993h . 

When available, 100 ^im flux densities from the IRAS FSC 
were also obtained for each 6dFGS-FSC source. For the ~ 
20 per cent of objects which did not have a 100 ^m detection the 
100 nm flux density was assumed to be twice the 60 ^m flux den- 
sity based on the calculated average (logfgioo um/geo um)) = 0.3 
with a scatter of 0.2 from the bright galaxy sample of Soifer et al] 
(J989). 60 and 100 /im flux densities are then converted into the 
quantity FIR (in units of Wm^^) defined by 

FIR = 1.26 X 10"" (2.58S'6o + ^loo ^m) (1) 

where 5*60 jim and Sioo Mm are in lanskys faelou et al] Il985l) . 
FIR is a measure of the total far-infrared flux between 42.5 and 
122.5 fim. 
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4.2 The FIR-radio ratio q 

The radio-FIR correla tion is often para metrised by the FIR-radio 
flux ratio parameter q ( lHelouetalJl985h . This is defined by 



= log 



/ FIR/(3.75 X 10' 



Si.. 



(2) 



where FIR (defined by equation[T) is divided by the factor 3.75 x 
10^^ Hz (the frequency at 80 ^l.m) to convert FIR to Wm"^ Hz"^ 
The mean value of (q) = 2.28 with rms scatter Oq — 0.22 for all 
2411 objects with measured flux densities at both 60 and 100 /jm . 
aq is somewhat larger than the value found bv lCondon et alj j2002l) 
for a more nearby sample of UGC-NVSS galaxies. This discrepancy 
reflects the large uncertainties associated with fainter IRAS FSC flux 
densities as well as the larger radio-loud AGN population contained 
in the present sample. For the subset of 2 242 star-forming galax- 
ies with measured 60 and 100 /im flux densities (gs p) = 2.3 with 
rms scatter crgp — 0.18, in close agreement with ICondon et all 
I2OO2) for the UGC-NVSS sample. This value also agrees well 
wit h results from sampl es of stronger IRAS sources in spiral galax- 
ies jCondonetal.ll99lh . For the subset of 169 radio-loud AGN with 
measured 60 and 100 /^m flux densities {qagn} ~ 2.0 with rms 
scatter ctagn ~ 0.5. The larger scatter is expected for radio-loud 
AGN as these are not expected to be correlated as strongly and the 
smaller average q value is caused by the stronger flux densities of 
radio-loud AGN. 

A cutoff of q — 1.8 is often used as a diagnostic to dis- 
tinguish between s tar-forming galaxies and radio-loud AGN (eg. 
ICondon et al.ll2002l) . and it appears that the spectroscopic classifi- 
cation of radio sources used in this work agrees well with this diag- 
nostic. Only 10 (< 1 per cent) of the SF galaxies in the combined 
6dFGS-NVSS-FSC sample have g < 1.8, the value below which 
a galaxy is 3 times more radio loud than the mean for star-forming 
galaxies. These 10 galaxies may come from a class of "compos- 
ite" radio source, whose radio emission is a mixture of both star- 
forming and radio-loud AGN activity or they may be candidate 
"taffy" galaxy pair s which recently suffered a direct collision (eg. 
ICondon et al.lll993l) . The SF spectral classification of these objects 
is retained in further analysis. 



4.3 Far-infrared colours 

Another diagnostic which can be used to distinguish be- 
tween radio-loud AGN and star-forming galaxies is the far- 
infrared spectral index between 60 and 25 ^im, qir = 
log(S'25/jm/5'60Mm)/log(60/25). ciiR is 3 measure of the temper- 
ature of the FIR emitting dust. Star-forming galaxies have cooler 
(T < 50 K) dust emission because of their extended star and gas 
distributions and therefore have values of am < —1.5. The central 
engines of AGNs heat the dust to warmer temperature s (T > 100 K) 
and have values of aia > — 1.5 l ldeGriipetalll992h . 

Sixty radio-loud AGN and 769 SF galaxies in the sample were 
detected in the IRAS FSC at 25 ^im. Fig.|9]plots am vs. q for these 
galaxies, comparing the dust temperature of each galaxy with its 
degree of radio excess. The vast bulk of SF galaxies lie in the lower 
right of the plot as expected from the radio-FIR correlation. The 
majority of radio-loud AGN have am > —1.5 (many lie outside 
the boundary on the upper left of Fig. |9]l indicating that the FIR 
spectral index, used in concert with the radio-FIR correlation agrees 
well with spectral classification as a discriminator of star-forming 
galaxies and radio-loud AGN. 
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Figure 9. apm vs. q (defined by equation [3 for the 769 star-forming 
galaxies (dots) and the 60 radio-loud AGN (circles) in the combined 6dFGS- 
NVSS-FSC sample which also had a detection in the IRAS-FSC at 25 /im. 
A dashed horizontal line is shown at apm = —1.5 to delineate the regions 
occupied by "warm" («fir > —1.5) and "cool" (apiR < —1.5) IRAS 
galaxies. The canonical value of g = 1. 8, used to distinguis h star-forming 
galaxies from AGNs in other surveys (eg. ICondon et alj2002h is shown as a 
dotted veilical line. 



4.4 Composite optical spectra and IRAS comparisons 

For about 90 per cent of the radio-detected 6dFGS galaxies in Ta- 
ble |3] a classification as either AGN or SF (as discussed in Sec- 
tion l2.5t appears straightforward. Most of the galaxies we have clas- 
sified as SF are detected by IRAS, and almost all of them lie in the 
lower right-hand quadrant of Fig. |9l as expected if their FIR emis- 
sion arises from dust heated in star-forming regions. In contrast, the 
IRAS detection rate of the objects classified as Aa in Table[3]is close 
to zero, as expected if their inner regions are largely free of gas and 
dust and their radio emission comes exclusively from an AGN. 

For the 10 per cent of galaxies in Table [3] which show op- 
tical emission lines with ratios characteristic of an AGN, and 
which are therefore classi fied as Aa e or Ae , the situation is less 
clearcut. As discussed by ISest et ah I l l2005bh . a possible problem 
with spectral classification of emission line AGN as radio-loud is 
that em ission-line AG N activity is often be accompanied by star for- 
mation jKauffmann e t al. 2003). Even if the AGN is radio-quiet, this 
associated star formation will give rise to radio emission. For any 
such galaxy included in the sample, the origin of the radio emission 
will be star formation, but the optical spectrum could still be dom- 
inated by emission lines from a (radio-quiet) AGN, leading us to 
classify the object as Aae or Ae. 

Although the present data are insufficient to disentangle the 
separate contributions from star formation and the AGN to the radio 
emission from these galaxies, we can estimate what fraction of them 
might have been misclassified spectroscopically by comparing our 
spectroscopic classification with the radio-FIR diagnostic. High res- 
olution Very Long Baseline Interferometry (VLBI) imaging of these 
sources may be used to distinguish the parsec-scale AGN compo- 
nents of these galaxies f rom the kpc scale radio structure resulting 
from star formation (eg. Norris e/a/.l ll992l) . Such observations are 
outside the scope of the present paper. 

Twenty spectroscopically classified AGN in the 6dFGS-NVSS 
sample lie in the lower right region of Fig.|9]and would be classified 
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Figure 10. The radio-FIR correlation of all 2 690 star-forming galaxies 
(dots) and 208 radio-loud AGN (open circles) in the combined 6dFGS- 
NVSS-FSC sample. A best fit line (equation[3]in the text) based on a subsam- 
ple of 2 242 star-forming galaxies with accurate 60 and 100 /im flux densities 
is also shown. 



as star-forming galaxies by the radio-FIR diagnostic. This imphes 
that about 1/3 of all spectroscopically classified emission-line AGN 
(Ae and Aae class) in the 6dFGS-NVSS sample (~ 250) fail the 
radio-FIR diagnostic. It is probable that some of these may be gen- 
uinely misclassified as radio-loud AGN from their optical spectra 
and derive the majority of their radio emission from star formation. 
Assuming that all absorption-line radio-loud AGN (Aa class) have 
been correctly classified, the disagreement between spectroscopic 
classification and the radio-FIR diagnostic for 6dFGS-NVSS AGN 
is ~ 10 per cent. We estimated a similar reliability for our classi- 
fication scheme in Section 12.51 from repeated visual inspection of 
the spectra. It is clear that difficulties in classifying emmision-line 
spectra with composite or borderline AGN-l-SF properties are the 
primary cause of misclassification of radio sources in our sample. 
We believe that this is a worst-case estimate of the reliability of our 
spectroscopic classification and note that on the whole the spectro- 
scopic classification technique agrees well with the radio-FIR cor- 
relation diagnostic. 

Only 38 per cent of objects in the 6dFGS-NVSS sample were 
detected in the 6dFGS-FSC sample. This makes the radio-FIR cor- 
relation diagnostic unsuitable for the sample presented in this pa- 
per because a non-detection of a galaxy in the 6dFGS-FSC sample 
does not preclude it from being a star-forming galaxy in the 6dFGS- 
NVSS sample. Deeper wide-angle far-infrared surveys are required 
to detect star-forming galaxies out to the redshifts probed by the 
current generation of optical redshift surveys. 



sample of ISadler et al.l j2002h . The FIR luminosities show less cor- 
relation with radio power for radio-loud AGN; these are found 
across the whole range of FIR luminosity sampled. Many of them 
have "excess" radio powers, causing them to lie above the correla- 
tion defined by the star-forming galaxies. These radio-excess IRAS 
galaxies are of interest for understa nding the connec tion between 
star formation and AGN in galaxies dOrake et al.ll2003h . 

A least-squares line of best fit to the subsample of 2 242 star- 
forming galaxies for which accurate 60 and 100 /im flux densities 
are available in the IRAS-FSC is also shown in Fig.[TO] The line has 
the form 



log(Pi., 



(1.06 ± 0.01) log(LFiR) + (11.11 ± 0.1) 



(3) 



which is similar to that found in other de rivations of th e radio - 
FIR correlation for brighter galaxies (eg. Condon et alj ( Il99lh ; 
ICondon & Broderickl ( ll988h ; lDevereux & Bales! ( Il989h ). 



4.6 Ultraluminous infrared galaxies 

Ultra-luminous Infrared Galaxies (ULIRGs) with Lfir > 10^^ Lq 
which are powered by a mixt ure of starburst activity and an 
AGN jSanders & Mirabell 1 19961) are of interest for understand- 
ing the starburst-AGN connection. There are two ULIRGs in the 
present sample; such objects are extremely rare in the local uni- 
verse. The nearest known (Arp 220) has redshift 0.018. The local 
U LIRG space density, e stimated from the FIR luminosity function 
of ISaunders et all l ll990l) . is < 10"** Mpc"^. The enclosed volume 
of the 6dFGS-NVSS sample (see Tabled is ~ 3.9 x 10* Mpc^ 
meaning that less than 4 ULIRGs a r e expe cted in the 6dFGS-NVSS 
volume. Interestingly, I Sadler et al.l ( l2002h found 7 ULIRGs in the 
2dFGRS-NVSS sample, in a volume of 5.0 x lO'' Mpc^ in the red- 
shift range 0.15 < z < 0.3, which implies their space density is 
1.4 X lO"'^ Mpc^^ within the more distant 2dFGRS redshift range. 
This suggests that these objects evolve strongly with redshift and 
is consiste nt with their measured evolution of (1 -|- z)^ between 
< 2 < 1 dCowie et alj2004l;ISandersll2004h . 



5 THE LOCAL RADIO LUMINOSITY FUNCTION 

The luminosity function ($(P)) of a sample of astronomical objects 
is a measure of the variation in their space density with luminos- 
ity. The local radio luminosity function (RLF) is t he global average 
space density of radi o sources at the present epoch jAuriemma et 
ll977l ; ICondonll989l) . An accurate derivation of the local radio lumi- 
nosity function is important as a present epoch benchmark for stud- 
ies of the population of radio sources at higher redshift so that their 
cosmological evolution can be determined. In this section we de- 
rive the local radio luminosity function (RLF) for the 6dFGS-NVSS 
sample as a whole and for the star-forming galaxy and radio-loud 
AGN subgroups separately. 



4.5 The FIR-Radio correlation of 6dFGS galaxies 

Fig. 1101 shows the radio-FIR correlation (Pi. 4 vs. I/pir) of 6dFGS- 
NVSS galaxies with a 60 /xm detection in the IRAS FSC. Spectro- 
scopically classified AGN are plotted as circles and SF galaxies are 
plotted as dots. The SF galaxies show a strong correlation between 
radio power and FIR luminosity extending over 4 orders of magni- 
tude though the scatter tends to increase above Lfir — 10^^ Lq 
as was seen for IRAS-FSC detected galaxies in the 2dFGRS-NVSS 



5.1 Completeness 

The RLF can be calculated directly but requires a sample which is 
complete to the limits of all the surveys from which it is derived. 
We therefore must define a subsample of the 6dFGS-NVSS sample 
which is complete to both its if-band magnitude limit and 1.4 GHz 
flux density limit. In the A'-band, the 2MASS XSC from which the 
6dFGS primary sample is selected is complete a nd reliable well be- 
low the K = 12.75 mag. hmit of the 6dFGS jjarrett et alj|2000l ; 
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Ijones et ailllOOd) . The 6dFGS primary sample may also miss pop- 
ulations of radio sources in the local universe either because their 
hosts appear stellar on optical plates (eg. QSOs, BL Lacs or com- 
pact galaxies) or because they are too blue in colour to be detected 
in the 6dFGS primary sample (eg. Seyfert galaxies). A selection of 
such objects were observ ed as additional ta rgets during the 6dFGS 
under programme id. 125 Jjones et alj2004l) . A full analysis of them 
will be the subject of a future paper. 

The sharp drop in the source counts of NVSS radio sources 
below 2.8 mjy (Fig.|5j indicates that the NVSS catalogue is incom- 
plete below this limit. We applied a 1.4 GHz flux density limit of 
2.8 mJy to the 6dFGS-NVSS sample when calculating t he RLF. The 
same flux density limit was used by Sadler et al.' (2002') to measure 
the luminosity function of 2dFGRS-NVSS radio sources. There are 
6961 6dFGS-NVSS galaxies with Si.dGHz > 2.8 mJy, K < 12.75 
of which 6667 have 6dFGSDR2 measured redshift z > 0.003. Of 
these, 3 997 were assigned spectral class SF, 2 652 were assigned 
AGN and 18 were unclassified. These 18 unclassified spectra with 
measured redshifts were classified on the basis of their measured ra- 
dio power: those with Pi. 4 < 10^^ WHz"^ (9) were classed as SF 
and those with Pi. 4 > 10^^ WHz"^ (9) were classed as AGN. 

We corrected the measured luminosity function for incomplete 
sampling of the celestial sphere by the 6dFGSDR2 and for spec- 
troscopic incompleteness of the 6dFGS-NVSS sample. To account 
for incomplete sampling of the celestial sphere we normalised all 
volumes by the effective area of the 6dFGS-NVSS sample derived 
in Section [231 Of the 6961 6dFGS-NVSS objects meeting our se- 
lection criteria, 280 had spectra which were too poor in quality to 
determine a redshift and 14 were associated with galactic stars. The 
poor quality spectra arose for a myriad of reasons which were pri- 
marily instrumental (eg. broken fibres, misplaced buttons etc.) and 
are expected to be a random subset of the data. The spectroscopic 
incompleteness of 4 per cent causes the luminosity function to be un- 
derestimated, so values of $ (Pi. 4) have been increased by 4 per cent 
to compensate. 

5.2 Calculating the luminosity function 

We have measured the radio luminosity function using the 1/Vinax 
method of Schmidt ( 1968). Knax is the maximum volume in which 
a galaxy will satisfy all of the sample selection criteria, which in 
the case of the 6dFGS-NVSS sample are 5'i.4GHz > 2.8 mJy, K < 
12.75 mag. and 2 > 0.003. 

We have corrected the measured RLF for galaxy clustering at a 
distance s centred on our own Galaxy using 

^ = 1 + (4) 
p 3 — 7s V s / 

( |Peebles|[l980l) . as discussed bv 'Condo n et alj (I2OO2 I). Here pp/p 
is the expected overdensity near our own Galaxy, or the space den- 
sity (pp) of local galaxies divided by the average space density (p) 
of all galaxies, 73 is the slope and so the correlation scale length 
from a power-law fit, ^ (s) = (s/so) of the two-point correla- 
tion function in redshift space. The two-point correlation function 
of local radio s ources has be en derived from a subset of the 6dFGS- 
NVSS sample dMauchll 20061 ), and is adequate to describe clustering 
of radio sources centred on our own galaxy. Derivations of 7s and 
So for radio sources in the 6dFGS-NVSS sample yielded values of 
7s = 1.57 and So = 10.07. (s) has a power-law form for distances 
s < 30Mpc, and ^ (s) ~ for s > 30Mpc. Therefore to correct 
for the local overdensity, we have multiplied the volume within s 
by equation|4]to calculate Knax for s < 30Mpc. In practise, V^max 
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Figure 11. The local luminosity function at 1.4 GHz measured from all the 
radio sources in the 6dFGS-NVSS sample. The curve is the sum of the con- 
tributions to the luminosity function from fits of equations |5]and[6]to the SF 
and AGN data respectively. 
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Figure 12. The local luminosity function at 1.4 GHz derived separately 
for the radio-loud AGN (circles) and star-forming galaxies (crosses) in the 
6dFGS-NVSS sample. The two curves are the fits of equations[5]and|6]to the 
SF and AGN data respectively. 



is always much larger than the local clustering volume, and so this 
correction has made little difference to the results. 



5.3 Results 

Table |5] lists the measured local radio luminosity function for the 
star-forming galaxies, radio-loud AGN and the 6dFGS-NVSS sam- 
ple as a whole. The 4006 star-forming galaxies have (V/Vmax) ~ 
0.509 ± 0.005, the 2 661 radio-loud AGN have (F/Fmax} = 
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Table 5. Local RLFs at 1 .4 GHz for the radio-loud AGN, star-forming galaxies and for all radio sources. 



All galaxies 



loglO -Pi. 4 
(W Hz-i) 



N 



log <!> 
(mag~^ Mpc~ 



Star-forming galaxies 
1V log $ 

(mag~l Mpc~3) 



Radio-loud AGN 



N 



log $ 
(mag^^ Mpc^-^) 



20.0 
20.4 
20.8 
21.2 
21.6 
22.0 
22.4 
22.8 
23.2 
23.6 
24.0 
24.4 
24.8 
25.2 
25.6 
26.0 
26.4 
Total 

(V/Vm 



3 
46 
116 
319 
654 
1266 
1496 
1138 
658 
378 
259 
183 
81 
49 
16 
3 
2 

6667 



-2.90 
-2.52 
-2.84 
-2.82 
-3.00 
-3.25 
-3.68 
-4.23 
-4.78 
-5.06 
-5.36 
-5.57 
-6.03 
-6.33 
-6.74 
-7.30 
-8.12 
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+0 
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+0 


27 


-0 


85 



3 
43 
103 
296 
589 
1106 
1119 
588 
133 
26 



'2.90lg:^° 
-2.56in? 

_2 SQ+0-04 



-2.85 
-3.05 
-3.31 
-3.79 
-4.45 
-5.36 
-6.12 



0.05 
+0.03 
0.03 
+0.02 
0.02 
+0.01 
0.01 
+0.02 
0.02 
+0.02 
0.02 
+0.04 
0.05 
+0.11 
0.14 



4006 
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+0 


27 


-0 


85 



0.518 ± 0.004 



0.509 ± 0.005 



0.532 ± 0.006 



0.532 ± 0.006 and all the 6 667 radio sources in the combined sam- 
ple have (V/Vmax) = 0.518 ± 0.004. The {V/Vra^^) value for 
radio-loud AGN is more than 3cr from the value of 0.5 expected 
if there were no significant clustering or evolution in the sample. 
Some evolution of the radio-loud AGN population is probable over 
the 1-2 Gyr lookback time of the 6dFGS-NVSS sample. The local 
radio luminosity function of all 6dFGS-NVSS galaxies is shown in 
Fig. 111! its statistical errors are of order 1 percent or less over 5 
decades of radio luminosity. Separate local radio luminosity func- 
tions for both star-forming galaxies and radio-loud AGN are shown 
in Fig. [T2] These cross over at Pi. 4 = 10^'^ WHz~^, star-forming 
galaxies dominate the population of radio sources below this power 
and radio-loud AGN dominate the population above it. 

Galaxy lumin osity functions are commonly fitted by the 
Schechter function jSchechteilll976h . This function turns over more 
steeply toward high luminosities than radio luminosity functions of 
star-forming galaxies. Instead the radio luminosity function is com- 
monly fitted by the parametric form given by 



HP) 



exp 



log.o (1 + P/P, 



(5) 



as was used to fit th e luminosity function of IRAS galaxies by 
ISaunders et alj ( Il990h . The best-fitting parameters of equation[5]for 
6dFGS-NVSS star-forming galaxies are: 



C 
P, 



10 



-2.83±0.05 



mag. ^ Mpc 



^q21.18±0.22^jj^-1. 

a = 1.02 ±0.15; 
cr = 0.60 ± 0.04. 

Equation [5] is inadequate to describe the RLF of radio-loud AGN 
and we have instead fitted these data with a 2 power-law function 
analogous to the optical luminosity function of quasars: 

C 



HP) 



[p./PY + {P*/P)^ 



(6) 



terown et all 1200 ll ; iDunlop & Peacocklll990h . The best fitting pa- 
rameters of equation |6] for 6dFGS-NVSS radio-loud AGN are: 

C = 10-^-^°±°-^^ mag.-^ Mpc-^ 

P. = 10'*■^^±''■='°WHz-^ 

Q = 1.27 ±0.18; 

P = 0.49 ± 0.04. 

The curves plotted in Fig. [12] show the fits of equations |5] and [6] 
to the star-forming galaxies and radio-loud AGN respectively. The 
radio luminosity function of all 6dFGS-NVSS galaxies is the sum 
of the two contributions from star-forming galaxies and radio-loud 
AGN, as shown in Fig.[TT] The curve in this figure is the sum of the 
two separate fits to the radio luminosity function for the star-forming 
galaxies and radio-loud AGN. 



5.4 Discussion 

The radio luminosity function of active galaxies maintains a power- 
law form for about 5 orders of magnitude before turning over above 
Pi. 4 > 10^*'^WHz~^. This power-law form of the radio-loud 
AGN luminosity function can only continue to fainter powers if it 
does not exceed the space density of the luminous galaxies in which 
radio-loud AGN preferentially reside. Taking the best fit Schechter 
function to the B-band luminosity fo r early-type galaxies in the 2dF- 
GRS from iMadgwick et"ai] j2002h . and assuming that radio-loud 
AGN typically live in the brightest early type galaxies (see Fig. [S] 
top), the space density of early-type galaxies with {Lb > Lg) is 
(9.3±0.3)xlO"'*Mpc-^ (forffo = 70 kms'^Mpc"^). The space 
density of radio-loud AGN with Pi. 4 > 10^°-^ WHz"^ from the 
6dFGS-NVSS luminosity function of Fig.[T2]is ?a 8 x 10"" Mpc'^ 
If the power-law form for radio-loud AGN continues to fainter radio 
powers a maximum must be reached at Pi. 4 ~ 10^^^ W Hz~^, just 
beyond the limit of the present data. Though there is some evidence 
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(eg. from Fig. [8} that radio-loud AGN of lower radio power may live 
in less luminous host galaxies which have higher space density. 

In recent models, in which the radio-loud AGN luminosity 
function is interpreted as the distribution of the time spent by an 
AGN at a given radio power, the strong variability of radio-loud 
AGN is rela ted to the heating of cooling flows in galaxies and 
clusters (eg. iNipoti & Binnevll20o3 : ISest et al]|2006l) . From a fit 
to the rad io luminosity function of 2dFGRS-FIRST radio sources, 
iNipoti ife Binnev (2005) predict that the AGN luminosity function 
will turn over at P1.4 = iQ^O-^io i WHz"\ just at the faint hmit 
of the 6dFGS-NVSS luminosity function. Such models are poorly 
constrained by observational data at present and require further con- 
straints on the relationship between the mechanical luminosity of 
heating of cooling flows by radio jets in AGN as well as a more 
detailed understanding of the X-ray luminosity function of AGN. 

5.5 Comparison of recent RLF measurements 

Fig. [13] shows a comparison of the RLFs of 6dFGS-NVSS galax- 
ies with those of UGC-NVSS galaxies and 2dFGRS-NVSS galax- 
ies. Values of the 2dFGRS-NVSS luminosity function have been 
shifted from their published values to Hq — 70kms^^ Mpc~^. 
The properties of these three surveys are compared in Table [T] 
We make no comparison of our measured RLF with that measured 
from SDSS-NVSS/FIRST galaxies because the normalisation of the 
SDSS-NVSS/FIRST RLF is tied to the 2dFGRS-NVSS RLF which 
it foll ows very closely o ver the complete range of radio powers sam- 
pled teest et al.ll2005bh . We also do not make any comparison with 
the 2dFGRS-FIRST RLF because the FIRST survey is insensitive 
to the large-scale radio e mission seen in many nearby radio sources. 
We therefore believe that lMagliocchetti et al.l ( l2002h have underes- 
timated the total flux density of many of the radio sources in the 
2dFGRS-FIRST sample. The 6dFGS-NVSS sample presented here 
samples a larger volume of space than the other radio source samples 
and probes lookback times intermediate between the nearby UGC- 
NVSS sample and the more distant 2dFGRS-NVSS one. These three 
luminosity functions are all drawn from the same radio source cat- 
alogue and differ only in the optical/near-infrared samples from 
which they are selected. The agreement between them is striking. 

5.5.1 Star-forming galaxies 

The 6dFGS-NVSS luminosity function of star-forming galaxies in 
the upper panel of Fig. [13] agrees with the other two over the 
range lO^"'^ < P1.4 < lO^^WHz"^ The falloff at the Pi. 4 = 
10^" WHz~^ point is caused by the small numbers of low lumi- 
nosity star-forming galaxies in the 6dFGS sample volume and can 
also be seen at the lowest radio powers in the luminosity function 
of 2dFGRS-NVSS galaxies. The UGC-NVSS luminosity function 
extends about an order of magnitude deeper in radio power for star- 
forming galaxies, reflecting the larger number of low-luminosity 
radio sources detected in a flux density limited sample covering a 
larger solid angle. The 6dFGS-NVSS luminosity function of star- 
forming galaxies does not fall as steeply as the other two at the 
highest radio powers (Pi. 4 > 10^^'^ WHz~^). We believe this is 
caused by the near-infrared selection of the 6dFGS-NVSS sample. 

5.5.2 Radio-loud AGN 

The radio-loud AGN luminosity function of 6dFGS-NVSS galax- 
ies agrees well with the UGC-NVSS radio-loud AGN luminosity 
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Figure 13. The local radio luminosity function at 1 .4 GHz for star-forming 
galaxies (top) and radio-loud AGN (bottom), comparing results from three 
differe nt samples. The UGC-NVSS luminosity functions of Condon et alj 
are s hown as stars. The 2dFGRS-NVSS luminosity functions of 
ISadler et alj i2002l) are shown as boxes and have been modified from 
Hq = 50kms~l Mpc~^ to Hq = 70kms^^ Mpc~^ by the scaling factor 
Pi. 4 (X and <I> oc H'^. The 6dFGS-NVSS luminosity functions are 

shown as circles with a line linking the points to highlight their position. 

function. The UGC-NVSS radio-loud AGN sample consists of only 
294 objects and is therefore not as well constrained. Most of the 
large error bars in the UGC-NVSS radio-loud AGN luminosity func- 
tion overlap with the 6dFGS-NVSS points. A power-law fit to the 
6dFGS-NVSS luminosity function below logiJPi.4 (WHz^^)] = 
25.0 is 

logio[$(Pi.4)] = (-0.61 ± 0.02) log[Pi.4] + (9.33 ± 0.53) (7) 

and to the 2dFGRS-NVSS luminosity function is 

logio[<l>(Pi.4)] = (-0.62 ± 0.03) log[Pi.4] + (9.42 ±0.70). (8) 

The 2dFGRS-NVSS luminosity function has the same power-law 
form as the 6dFGS-NVSS luminosity function over nearly four or- 
ders of magnitude. The present data are insufficient to measure any 
possible evolution over the narrow redshift range of 5 = 0.07 to 
5 = 0.14 between the two samples. 
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Figure 14. The variation in star-formation density (psf(SFR)) with star- 
formation rate for galaxies with star-formation rates between 0.1 and 
100 Mq yr~^ as derived from equation 1111 Values of $ have been taken 
from Table [5] for star-forming galaxies with the errors scaled appropriately. 
The curve is derived from the fit of equation |5]to the luminosity function of 
star-forming galaxies. 



6 THE STAR-FORMATION DENSITY AT THE PRESENT 
EPOCH 

Radio emission from star-forming galaxies is the product of mas- 
sive {M > 8AIq) stars whicli liave short lifetimes (r ~ 10 yr). 
Roughly 90 per cent of this radio emission comes from non-thermal 
synchrotron electrons, which have been accelerated in the remnants 
of type II supemovae, and the remainder comes from ionised hydro- 
gen gas in Hll regions heated by the most massive stars. This means 
that the radio luminosity of star-forming galaxies is roug hly pro- 
portio nal to the rate of recent (r ^ 10* yr) star-formation jCondonI 
Il992h . It is therefore possible to use the star-forming galaxies in this 
sample to estimate the present epoch star-formation density of the 
universe for galaxies with K < 12.7 5 jnag., which is the zero point 
of the Madau diagram ( Mada u et al. Hl996h . Estimates of the local 
star-formation density made from radio surveys are particularly ro- 
bust as they are free from the effects of dust extinction. Effects of 
extinction are minimised in the 6dFGS-NVSS sample, as the radio 
selection is made from a near-infrared galaxy sample. 

6.1 Calculating the star-formation density 

It is possible to estimate the local star-formation density as a func- 
tion of star-formation rate directly from the radio lu minosity func- 
tion u sing a method described by Cram ( 1998) and iHaarsma et al] 
Assuming a typical Salpeter-like initial mass function of the 

form 

^(Af) oc M"^-^'^ (9) 

between 0.1 and lOOM©, the star-formation rate of stars more mas- 
sive than 0.1 M0 (in Mq yr^^) can be calculated from the 1.4 GHz 
radio power by the relation 

jSullivan et alj|200l]) . The star-formation density can then be deter- 
mined as a function of star-formation rate by multiplying the star- 



formation rate by the space density of radio sources, or simply 

Psf(SFR) = SFR(Pi.4) X <l>(Pi.4) (11) 

in Mq yr~^ mag.~^ Mpc~^. Fig. [14] shows equation [TT] for the 
luminosity function of 6dFGS-NVSS star-forming galaxies. Star- 
forming galaxies in the 6dFGS-NVSS sample trace star-formation 
rates between 0.1 and 500 Mq yr~^, though the major contribution 
to the local star-formation density comes from galaxies with star- 
formation rates around 10 Mq yr~^ as has been seen in other de- 
termi nations of the star-f ormation density of the local u niverse (eg. 
radio:'Sa dler et alj l l2002h and Hai lGallego et al.l ( Il995h ). 

Integrating underneath the curve in Fig. [14] gives an estimate 
of the global star-formation density at the present epoch (psp)- The 
vast majority of radio sources in this sample contribute to the inte- 
gral near its maximum at 10 Mq yr~^ ; the less certain points below 
this maximum make little contribution to the integral. This implies 
that the global star-formation density will be accurately constrained 
by the data. Rather than compute the global star-formation density 
from the binned data in Fig. [14] it is preferable to compute psF di- 
rectly from a sum over each sample galaxy. This is done by comput- 
ing the local power density function (u(Pi,4)), 

The total radio power produced per unit volume of space Usf 
(in units of WHz"'^ Mpc""^) is just the integral over all radio powers 
of the power density function 

Usf = / u{PiA)dPiA, (12) 
Jo 

which can be calculated directly as the sum over all galaxies of 
Pi.4/Vmax. Usf is directly proportional to the star-formation rate 
density psp and is computed using the relation 

PSF = 1.13 X 10-^'[/sF (13) 

the factor 1.13 x 10^^^ comes from the conversion between star- 
formation rate and radio power in equation [Tol 



6.2 Results 

The global power density of star formation at the present epoch 
implied by the 6dFGS-NVSS sample is ;7sp = (1.91 ± 0.09) x 
10^'^ WHz~^ Mpc~^. This translates to a global star-formation den- 
sity 

PSF = (0.022 ± 0.001) Mq yr"^ Mpc"'' (14) 

The uncertainty in psF quoted here is purely statistical and takes no 
account of the errors arising from the conversion from Pi. 4 to star- 
formation rate in equation [TO] Errors in this conversion come pri- 
marily from predicting the star-formation rate from the Type II SNe 
rate, which contribu tes 90 per cent of the radio power at 1.4 GHz 
( ISullivan et al.l20()II) . from models of the total radio spectral energy 
per Type II supernova and from the assumption of a Salpeter-like 
initial mass function in equation|9] 

Values for psF measured from the UGC-NVSS and 2dFGRS- 
NVSS samples are shown in Table [T] The value of psF for 
6dFGS-NVSS gal axies is intermediate between values those for 
the other samples. ICondon et ai] l l2002h derived psF ~ 0.018 ± 
0.001 Mpyr-^Mpc "^ for the UGC-NV SS sample of 1672 star- 
forming galaxies and lSadler et al.l ( I2OO2I) derived psF = 0.031 ± 
0.006 Mq yr"^ Mpc"^ (converted to 77o = 70 km s"^ Mpc"^) for 
242 star-forming galaxies from the 2dFGRS-NVSS sample. Dif- 
ferences in the measured values of psF between the three radio 
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samples are likely to be caused by their different optical/near- 
infrared magnitude limits. Fainter and hence less massive galax- 
ies can make a significant contribution to the local star formation 
density ijB rinchmann et al. 2004). For example, from Fig. 15 of 
iBrinchmann et alj (|2004) we estimate that the K — 12.75 mag- 
nitude limit of the 6dFGS-NVSS means the sample misses at least 
~ 0.004 Mq yr~^ Mpc"'' of the local star-formation density mea- 
sured from Ha. The 2dFGRS-NVSS sample on the other hand has 
an optical limit about 2 magnitudes fainter than the 6dFGS-NVSS 
sample and misses only a negligible amount of local star forma- 
tion. Further differences between the three samples are probably the 
result of evolution in the star-formation density over the narrow red- 
shift range. 

Uncertainties in the relationship between monochromatic radio 
power and star-formation rate do not affect models describing the 
relationship between star-formation rate and luminosities at shorter 
wavelengths. However, dust extinction within the Hll regions where 
star formation occurs contributes increasingly to the model uncer- 
tainties at these shorter wavelengths. The value of psF measured 
here has been compared with other rec ent determinations measure d 
from a broad range of wavelengths bv ' Hopkins & BeacomI ( l2006h . 
The interested reader is referred to this paper for more information. 
The general agreement (usually to within a factor of 2) between val- 
ues of psF at different wavebands is reassuring. 



7 BIVARIATE RADIO-NIR LUMINOSITY FUNCTIONS 

The local RLF of radio-loud AGN presented in section |5] contains 
large numbers of galaxies in most bins of radio power, which en- 
ables us to split it into statistically significant subclasses. In partic- 
ular, the local radio luminosity functions can be broken into bins 
of absolute K-band magnitude (Mk) to examine the radio prop- 
erties of galaxies of different near-infrared luminosities. The near- 
iirfrared selection of the 6dFGS-NVSS sample is of particular im- 
portance here, as the iC-band galaxy luminosity is an extinction- 
free estimator of the mass of the old stellar population in galaxies. 
The vast majority of radio-loud AGN reside in early type galax- 
ies whose absolute K-band magnitude gives an indication of their 
bulge mass which in tu rn is correlated with the black hole mass (eg. 
|m arconi & Hunt 2003). It is therefore possible to examine how the 
properties of radio-loud AGN may vary with black-hole mass. 

One of the fir st bivariate radio - optica l luminosity functions 
was determined by lAuriemma et alj ( 1 19771) from a sample of ra- 
dio detections of bright elliptical galaxies. They showed that the 
probability of an elliptical galaxy to host a radio-loud AGN in- 
creases strongly as a function of optical luminosity. Further stud- 
ies using larger sa mpl es of elliptical ga l axies were carried out by 
ISadler et al.' (1989') and'Ledlo w & Owed ( Il996l) who verified the re- 
sult of Auriemma et al. ( 1977) and found that the break in radio lu- 
minosit y function s caled strongly with host galaxy luminosity. More 
recentlv lBest et alj fcoOSa) used the SDSS-NVSS/FIRST sample of 
2 215 galaxies to calculate accurate bivariate radio luminosity-mass 
functions, from both host-galaxy mass and black-hole mass. They 
found that the fraction of galaxies hosting a radio-loud AGN in- 
creases strongly a s a function of both galaxy and black hole masses. 
iBest et alj J2005al) also argued that there was no evidence for any 
dependence of the break luminosity on black hole mass. 

Here, we measure the bivariate radio-NIR luminosity function 
for the first time which we can compare with bivariate radio-optical 
luminosity functions. The 6dFGS-NVSS bivariate luminosity func- 
tion is calculated from accurate observed quantities measured from 
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Figure 15. The fractional luminosity function of radio-loud AGN calculated 
using the procedure described in section lTTl of the text. Fk{Pia) is calcu- 
lated in 4 bins of absolute K magnitude (Mx ) as labelled in the figure. The 
two Mx bins centred at —25.5 and —26.5 are brighter than = —24.6. 

homogeneous all-sky surveys. We therefore avoid any errors which 
may arise from calculating the bivariate luminosity function from 
derived quantities. For examp le the large observe d scatter in the 
A/bii — o" relation of 0.3 dex tTremaine et alj|2002l) is a significant 
source of error in the biva riate radio l uminos ity versus black-hole 
mass function measured bv lBest et al.U2005al) . 

7.1 Calculating the bivariate luminosity function 

When separating the radio luminosity function into Mk bins it is 
preferable to define a normalised, or fractional luminosity function 

pK 

which is the radio luminosity function of galaxies in Mk bin 
Mk ± AMk (0if(Pi.4) in Mpc"^) divided by the volume den- 
sity ipK in Mpc~^) of all galaxies with absolute magnitude in the 
Mif -magnitude bin Mk ± AA//f . The radio luminosity function 
in each _ftr-band magnitude bin 4>k{Pia) is calculated using the 
1/Knax method ( ISchmidl|[l968l) in 1 -magnitude radio luminosity 
bins of width 0.4 in Iog(Pi.4). This gives 4>k{,Pia) in units of 
mag~^ Mpc~'^. With this definition FK(Pi.4,)dP\A represents the 
fraction of all galaxies in the absolute near-infrared magnitude bin 
Mk ± A Mk that have radio power bet ween Pi . 4 and Pi . 4 -I- dPi . 4 . 

To determine pK, the volume density of galaxies with absolute 
near-infrared magnitude Mk, it is convenient to use the K-bwaA 
galaxy luminosity function It is important that the A"-band lu- 
minosity function used be subject to the same selection criteria as 
the radio sample, so as not to introduce biases from cosmic variance 
into the data. We use the 1/VInax det ermination of the 6 dFGS K- 
band luminosity function described bv llones et alj 1 I2OO6I) . Its error 
bars are much smaller than those in the 6dFGS-NVSS luminosity 
function and are therefore ignored when calculating Fk{Pia)- 

7.2 The bivariate radio-NIR luminosity function of 
radio-loud AGN 

Table [6] lists and Fig. 1 151 shows Fk{Pia) for radio-loud AGN in 
four bins of Mk- Below Pk(Pi.4) = 0.1 the broad break between 
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Table 6. The bivariate radio-ii'-band luminosity function of radio-loud AGN. 
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seen in the RLF in Fig. [12] can 



be seen in the brightest Mk bins (Ma ~ —25.5, —26.5). In the 
faintest 2 bins (Mk ~ —23.5, —24.5), the data do not extend to 
high enough radio powers to detect a break. In the brightest bins of 
Mk the fractional luminosity function turns over below the radio 
power at which Fa-(Pi.4) = 0.1. It is probable that this turnover is 
caused by the normalisation of the fractional luminosity function; a 
larger fraction of radio sources cannot be measured at a given Mk 
because nearly all available galaxies are detected. The turnover oc- 
curs at fainter radio powers in less luminous galaxies, this is prob- 
ably because both less luminous galaxies and less powerful radio 
sources have higher space density. 

There is a clear decrease with radio power in the fraction of 
galaxies which are radio-loud AGN for all Mk, reflecting the de- 
crease in space density of radio-loud AGN with increasing radio 
power. For a given radio power, there is a strong increase in the 
fraction of galaxies which are radio-loud AGN with brightening 
Mk- More luminous galaxies are therefore more likely to host a 
radio-loud AGN as has been seen in previous determinations of the 
fractional radio-optical luminosity function of elliptical galaxies. To 
properly quantify this effect we determine the fractional luminosity 
function in its integral form. 
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Figure 16. The fractional luminosity function of radio-loud AGN in integral 
form. Fk(> Pia) denotes the fraction of all galaxies which host a radio- 
loud AGN with radio power greater than Pi. 4. Errors have been calculated 
using the bootstrap method with 10'' iterations. 



7.3 The integral bivariate radio-NIR luminosity function of 
radio-loud AGN 

The integral form of the fractional luminosity function at a particular 
radio power {Fk{> Pia)) is calculated in each Mk bin from the 
sum of Fk{Pia) between the largest bin of radio power (P™!"') to 
Pi . 4. It measures the fraction of all galaxies which are radio sources 
with powers greater than Pi. 4 in each Mk bin. It is not straightfor- 
ward to calculate the error at each point of Fk(> Pia) because the 
points are not independent. For this reason we calculated errors in 
each radio power bin of Fk{> Pia) using bootstrap analysis with 
lO*" iterations. 

Fig. [16] shows the fractional luminosity function of radio-loud 
AGN in integral form in four bins of Mk, with bootstrap errors. 
With the available data, radio powers at which 40-50 per cent of 



galaxies are detected as radio-loud AGN are reached at the brightest 
A'Ik bins. The function has approximately the same shape for all 
bins of A/a though this breaks down above Fk{> Pia) = 0.3 
which is caused by the tur nover in Fk{Pia) discussed in Sec- 
tion IT2I iBest et alj ( l2005j) also found a ceiling in the bivariate 
radio luminosity-mass function of SDSS-NVSS/FIRST galaxies at 
the 30 per cent level. We determine the dependence of /radio-ioud 
(the fraction of all galaxies that are radio loud) on La by scaling 
Fk{> Pia) to remove the dependence on Mk- The best-fitting 
scaling factor for Fk{> Pia) < 0.3 is /radioioud ex -La^; clearly, 
the probability of a galaxy to host radio-loud AGN is a strong func- 
tion of its _ftr-band luminosity. 

Given that the majority (> 70 per cent) of the radio-loud AGN 
in the 6dFGS-NVSS sample reside in early-type galaxies, their total 
A'-band luminosity should be close to or equal to their near-infrared 
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bulge luminosity (Lx bui)- Assuming this, the A"-band luminosi- 
ties of these radio-loud AGN are indicative of the black-hole mass 
of th e galaxy via the MBH-ijf.bui correlation arcom &Hunll 
l2003h . The precise scaling is uncertain though, as it is unclear at 
this stage of the precise contribution from disks to the Tf-band lumi- 
nosity of the radio-loud AGN in our sample. An accurate derivation 
of the if-band bulge luminosity of all 6dFGS galaxies is beyond 
the scope of this paper Accurate velocity dispersions will be mea- 
sured f rom spectra of the 15 000 brightest elliptical galaxies in the 
6dFGSl ljones et al.ll2004h . We expect to detect about 30 per cent of 
these in the NVSS (see Fig. |16t . A full analysis of these galaxies 
will be the subject of a future paper. 



8 CONCLUSION 

In this paper we have derived a catalogue of 7 824 radio sources 
from the NVSS catalogue which have been observed in the 
6dFGSDR2 if-selected primary sample. The main results of this 
paper are: 

• The median redshift of all radio sources in the sample is 
z = 0.046, which is similar to that of the 47 317 A'-band selected 
galaxies from which the radio sources were selected. 60 per cent 
of 6dFGS-NVSS primary targets are spectroscopically identified as 
star-forming galaxies which form a much more nearby population 
{z = 0.035) and 40 per cent of 6dFGS-NVSS primary targets are 
identified as radio-loud AGN which form a more distant population 
(2 = 0.073). 

• Star-forming galaxies and radio-loud AGN have quite distinct 
distributions in the plane of radio power vs. absolute K magni- 
tude. Star-forming galaxies tend to have radio powers Pi. 4 < 
lO'^^ WHz^^, though they have a wide range of absolute K magni- 
tudes. The radio-loud AGN span a wide range in radio power but are 
almost all found in the most luminous near-infrared galaxies with 
Mk > M*ji- 

• 2 690 star-forming galaxies and 208 radio-loud AGN were de- 
tected at 60 /im in the IRAS Faint Source Catalogue. All of the star- 
forming galaxies with IRAS-FSC detections in the 6dFGS-NVSS 
sample are found to lie on the radio-FIR correlation, with aver- 
age FIR-radio flux ratio parameter (qsf) = 2.3 with rms scatter 
asF = 0.18, in close agreeme nt with the results o f other recent 
radio-FIR surveys (UGC-NVSS: ICondon et ai]|2002l) . Many of the 
208 radio-loud AGN with IRAS-FSC 60 ^lm detections also lie on 
the radio-FIR correlation. However, these tend to have hotter IRAS 
60-25 fim colours. The combination of the radio-FIR correlation and 
the 60-25 pm colour as a means of AGN/Star-forming classification 
for radio sources agrees with spectroscopic classification for more 
than 90 per cent of all galaxies, which is reassuring considering the 
differences in these classification techniques. 

• We have measured accurate and homogeneous local radio lu- 
minosity functions at 1.4 GHz of star-forming galaxies and radio- 
loud AGN. The luminosity functions presented here agree well with 
recent determinations of the local radio luminosity functio n de- 
rived f rom different opti cal samples (eg. UG C-NVSS; Condon et al.l 
( l2002h . 2dFGRS-NVSS: ISadler et al] |20o2)). 

• The star-formation density at the present epoch has been mea- 
sured from the 6dFGS-NVSS sample and was found to psF = 
(0.022 ±0.001) Mq yr^^ Mpc^^, in agreement with previously de- 
rived values at radio and other wavelengths. 

• We have split the radio luminosity function of radio-loud AGN 
into four bins of Mk and compared with the K-hwA luminosity 



function of all 6dFGS galaxies to calculate a fractional luminosity 
function Fk{Pia)- We find that Fk{Pia) for radio-loud AGN in- 
creases strongly as a function of near-infrared luminosity, which in 
turn is indicative of an increasing probability for a galaxy to host a 
radio-loud AGN with increasing black-hole mass. 
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